UNCLASSIFIED 


AD  295  108 

Hep^duced. 
kf  Uie 

ARMED  SERVICES  TECHNICAL  INFORMAHON  AGENCY 
ARLINGTON  HALL  STATIffll 
ARLINGTffll  12.  VIRGINU 


NOTICE:  When  government  or  other  drawings,  specl> 
flcatlons  or  other  data  are  used  for  any  purpose 
other  than  In  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Govenunent  thereby  Incurs  no  responsibility,  nor  any 
obligation  tdiat soever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  In  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  Is  not  to  be  regarded  by  Implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rl^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  Invention  that  may  In  any  way  be  related 
thereto. 


295  108 


I 

I 

1 


%  ANALYTICAL  STUDY 

OF 

US  BEAM  CHARACTERISTICS 


•« 

V 


I 


by 

THE  BENDIX  CORPORATION 
ECLIPSE-PIONEER  DIVISION 

FA  A  Contract  No:  ARDS -451 
FA  A  Project  No:  114-1312D 
E-P  Project  No:  AD1091 


August  31,  1962 


FEDERAL  AVIATION  AGENCY  SYSTEMS  RESEARCH  AND  DEVELOPMENT  SERVICE 


SYSTEMS 


Analytical  Study 
of 

ILS  Beam  Characteristics 


by 

THE  BENDIX  CORPORATION 
ECLIPSE-PIONEER  DIVISION 


FAA  Contract  No:  ARDS- 451 
FAA  Project  No:  114-1312D 
E-P  Project  No:  AD1091 


August  31,  1962 


FEDERAL  AVIATION  AGENCY 
RESEARCH  AND  DEVELOPMENT  SERVICE 


NOTICE 


This  report  has  been  prepared  by  The  Bendix  Corporation's 
Eclipse-Pioneer  Division  for  the  Systems  Research  and  Development 
Service  (formerly  Aviation  Research  and  Development  Service), 
Federal  Aviation  Agency  under  Contract  Number  ARDS- 451,  Project 
Number  114-1312D.  The  contents  of  this  report  reflect  the  views  of 
the  contractor,  who  is  responsible  for  the  facts  and  the  accuracy  of 
the  data  presented  herein,  and  do  not  necessarily  reflect  the  official 
views  or  policy  of  the  FAA. 


THE  DENDIX  CORPORATION 
ECLIPSE-PIONEER  DIVISION 
TETERBORO.  NEW  JERSEY 


Analytical  Study 
of 

ILS  Beam  Characteristics 
for  the 

Systems  Research  and  Development  Service 
Federal  Aviation  Agency 


FAA  Contract  No:  ARDS- 451 
FAA  Project  No:  114-1312D 
E-P  Project  No:  AD-1091 


August  31,  1962 


Prepared  by: 


Project  Manager 


Reviewed  by: 


K.  Moses 
Senior  Engineer 


Approved  by: 

W.  A.  Platt 

Assistant  Chief  Engineer 


A.  Noxon 
Director  of  Advanced  Systems 
Development  Department 


FOREWORD 


This  report  was  prepared  by  the  Eclipse-Pioneer  Division 
of  The  Bendix  Corporation  as  required  by  the  Federal  Aviation 
Agency  under  Contract  FAA/ARDS-4.^1  ,  Project  Number  114-1312D. 

The  entire  study  was  administered  by  the  Advanced  Systems 
Development  Laboratory  with  J.  Doniger  as  Project  Manager.  Messrs. 
R.  Kostanty,  B.  Hanisch,  F.  Belsky,  and  A.  Flitt  performed  the  ana¬ 
log  computer  studies,  and  J.  Clair  studied  the  applicable  decision  mak¬ 
ing  techniques. 

The  cognizant  FAA  Project  Manager  was  Mr.  N.  J.  Proferes 
of  the  Aviation  Research  and  Development  Service;  who  together  with 
Mr.  D.  J.  Sheftel  and  Mr.  A.  B.  Winick,  provided  valuable  guidance 
throughout  the  project. 


The  Bendix  Corporation,  Eclipse-Pioneer  Division,  Teterboro,  N.  J. 


ANALYTICAL  STUDY  OF  ILS  BEAM  CHARACTERISTICS  FOR 
THE  SYSTEMS  RESEARCH  AND  DEVELOPMENT  SERVICE. 
August  31,  1962.  90  pages,  including  50  illustrations. 


ABSTRACT 


This  report  presents  the  results  of  analog  computer  studies 
and  analyses  of  the  noise  component  in  the  ILS  glide  path  and  localizer 
course  structures  due  to  terrain  reflections.  Actual  recordings  of  the 
noise  at  the  25  ILS  sites  which  will  probably  be  used  in  the  initial  lower 
minimum  operations  below  200  feet  were  furnished  to  the  contractor  by 
the  FAA.  The  studies  included  the  comparison  of  various  possible  ILS 
course/path  noise  figures  of  merit  as  they  related  to  the  performance  of 
presently  accepted  airborne  coupler  and  autopilot  systems  with  the  KC- 
135,  DC-7,  and  TF-102  aircraft. 

An  attempt  was  made  to  characterize  the  measured  ILS  course/path 
noise  in  terms  of  power  spectral  density  as  well  as  amplitude  and  to  use 
such  a  characterization  to  accept  or  reject  installations.  Only  weak  corre¬ 
lations  were  found  to  exist  between  aircraft  performance  and  the  possible 
power  spectral  figures  of  merit. 

A  satisfactory  course/path  noise  amplitude  specification  of  5  |xa 
in  the  localizer  and  20  pa  in  the  glide  path  at  low  altitudes  can  be  speci¬ 
fied  with  the  increased  probability  of  rejection  of  some  ILS  facilities  for 
approaches  below  the  present  200  feet  minimum. 

The  best  technical  solution  of  the  acceptability  criteria  problem 
was  found  to  be  the  application  of  the  measured  noise  to  a  complete  sim¬ 
ulation  of  several  models  of  airborne  systems. 
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SUMMARY 


This  is  the  compilation  of  two  final  reports  under  Phase  I  and 
Phase  lA  of  the  work  statement  specified  in  Contract  FAA/SRDS-451  as 
amended  July  1962. 

In  Phase  I  of  this  project  an  attempt  was  made,  using  actual  re¬ 
cordings  of  25  major  ILS  facilities,  to  find  correlation  between  test  air¬ 
craft  lateral  dispersions  and  selected  localizer  course  noise*  power  spec¬ 
tral  parameters  applied  to  such  facilities.  No  significant  correlation  was 
found  to  exist.  Only  the  RMS  aileron  and  roll  dispersions  were  found  to 
be  weakly  correlated  with  course/path  noise  power.  At  best,  these  trends 
may  provide  a  means  of  using  course/path  noise  power  spectra  in  a  lim¬ 
ited  fashion  as  a  secondary  criterion  for  localizer  course  acceptance  in 
addition  to  the  present  localizer  course  amplitude  criteria.  However,  the 
correlation  between  aircraft  performance  and  course  noise  power  improved 
with  the  use  of  mean  square  values  of  roll  and  pitch  attitude  or  lateral  and 
vertical  deviation  in  lieu  of  using  instantaneous  values  of  aircraft  position. 

Analysis  of  the  glide  path  data  revealed  that  the  path  noise  is  not 
statistically  stationary  in  nature  and  that  power  spectral  techniques  con¬ 
sequently  cannot  be  used. 

Approach  couplers  which  use  a  derived  rate  term  to  damp  the 
track  response  are  more  susceptible  to  the  adverse  effects  of  course/path 
noise  than  those  couplers  which  use  attitude  damping.  Heading  damped 
localizer  couplers  and  attitude  damped  glide  path  couplers  reduce  the  ef¬ 
fects  of  course/path  noise  by  almost  three  to  one. 

The  TF-102  lateral  responses  to  course/path  noise  were  two  to 
three  times  more  severe  than  the  KC-135  or  DC-7  responses.  This  is 
due  to  the  broader  band  airframe  response  as  a  result  of  the  low  moment 
of  inertia  about  the  roll  axis  and  the  higher  approach  speed  of  the  TF-102. 


'■‘The  term  "noise"  or  "course/path  noise"  used  in  this  report  only  re¬ 
lates  to  the  irregularities  in  the  ILS  course/path  structure  due  to  ter¬ 
rain  reflections. 
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The  TF-102  pitch  dispersions  did  not  significantly  differ  from 
the  other  aircraft  since  the  autopilot  configuration  masks  the  differ¬ 
ences  in  the  airframe  responses. 

The  employment  of  a  simple  filter  technique  to  accept  ILS  lo¬ 
calizer  and  glide  path  performance  did  not  result  in  sufficient  corre¬ 
lation  between  the  output  responses  of  relatively  simple  filters  and  the 
complete  ILS  system  simulation.  This  was  due  in  part  to  the  lower  or¬ 
der  approximation  of  the  filter  itself  to  the  actual  transfer  functions,  as 
well  as  to  the  non-Gaussian  and  nonstationary  characteristics  of  the  ac¬ 
tual  ILS  course/path  noise. 

An  arbitrary  course/path  noise  variance  technique  to  determine 
acceptability  of  ILS  facilities  did  provide  reasonable  correlation  be¬ 
tween  the  roll,  pitch,  vertical  and  lateral  deviations  using  a  random 
noise  generator  and  the  actual  course/path  noise  in  the  recordings. 

The  discrepancies  are  due  to  the  differences  between  the  peculiar  char¬ 
acteristics  of  the  actual  course/path  noise,  as  noted  above,  and  the 
characteristics  of  the  almost  purely  random  generator;  as  well  as  the 
fact  that  the  actual  noise  variance  does  not  change  with  range  in  the 
same  manner  as  the  generated  noise  was  changed. 

In  order  to  minimize  the  possibility  of  accepting  an  unsatisfac¬ 
tory  facility  for  low  approach  use,  the  results  of  the  course/path  noise 
variance  program  are  as  follows,  based  on  the  desire  to  restrict  95% 
of  the  aircraft  deviations,  due  to  this  noise,  in  the  region  of  the  middle 
marker  and  the  runway,  to  less  than  2  degrees  of  roll  and  pitch  attitude, 
lateral  displacements  to  less  than  15  feet  and  vertical  displacements  to 
less  than  4  feet; 

1)  The  localizer  course  noise  due  to  terrain  reflections  in  the 
region  between  the  middle  marker  and  the  runway  should  not 
exceed  5  microamperes  with  95%  probability.  *  The  course 
noise  in  the  region  of  the  outer  marker  can  be  30  microamperes 
(95%  probability  level)  without  deteriorating  performance  at  the 
lower  altitudes.  A  linear  relationship  between  the  allowable 


See  Section  3.  2.  3  for  description  of  the  test  procedure  to  be  used  in 
verifying  that  an  ILS  facility  is  within  the  limits  of  the  course/path 
noise  criteria. 
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course  noise  level  and  the  distance  from  the  landing  threshold  is 
recommended. 

2)  The  glide  path  noise,  due  to  terrain  irregularities,  from  the 
region  of  the  middle  marker  to  the  runway  should  not  exceed  20 
microamperes  with  95%  probability.  The  noise  in  the  region  of 
the  outer  marker  can  be  30  microamperes  f,95%  probability  level) 
without  deteriorating  performance  at  the  lower  altitudes;  with  a 
linear  relationship  between  the  allowable  path  noise  level  and  the 
distance  from  the  landing  threshold  being  recommended. 

The  foregoing  recommendations  represent  an  acceptable  solu¬ 
tion  since  they  can  be  implemented  relatively  simply  and  immediately. 
This  relatively  easy  solution  also  has  a  penalty  in  that  in  trying  to  pre¬ 
vent  the  possibility  of  falsely  accepting  an  unsatisfactory  facility  some 
acceptable  facilities  may  be  rejected.  This  study  has  shown  that  the  use 
of  increasingly  weighted  or  averaged  noise  characteristics  can  only  in¬ 
crease  the  uncertainty  between  the  predicted  and  actual  performance. 

The  best  technical  solution  short  of  actually  flying  a  series  of 
aircraft  and  automatic  systems  that  represent  the  types  that  will  use 
the  ILS  facility,  to  determine  whether  the  performance  is  acceptable,  is 
to  apply  a  recorded  version  of  the  ILS  course/path  noise  corrected  to  ac¬ 
count  for  the  motion  of  the  test  aircraft  to  complete  analog  computer  sim¬ 
ulations  of  the  various  systems  that  will  use  the  facility.  This  is  the  only 
technique  that  can  result  in  one-to-one  correspondence  between. predicted 
and  actual  performance. 
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1.  INTRODUCTION 


This  document  is  a  compilation  of  the  final  reports  of  Phase  I 
and  Phase  lA  of  the  ILS  localizer  course  and  glide  path  study  program 
outlined  inFAA/SRDS-451.  The  primary  objective  of  the  study  was  to 
determine  a  figure  of  merit  for  ILS  course/path  noise  which  provides 
a  strong  correlation  with  the  acceptability  of  the  approach  process.  In 
this  report  the  term  "beam  noise"  or  course/path  noise"  only  relates 
to  the  irregularities  in  the  localizer  course  and  glide  path  due  to  terrain 
reflections  of  the  carrier  and/or  sideband  energy  of  the  transmitted  sig¬ 
nals.  The  vector  sum  of  the  direct  and  reflected  course/path  signal  com¬ 
ponents  results  in  a  localizer  or  glide  path  reference  line  which  is  not 
straight  but  which  consists  of  irregular  granularities.  As  the  aircraft, 
which  is  involved  in  the  approach,  moves  through  these  regions  of  signal 
variations  along  the  course  line,  the  irregularities  in  the  course/path 
structure  cause  the  receiver  output  signal  to  fluctuate.  The  fluctuations, 
or  course/path  noise,  observed  in  the  aircraft  therefore  depend  on  the 
spacial  distribution  of  the  beam  irregularities  and  the  speed  with  which  the 
aircraft  traverses  them.  Automatic  ILS  coupler  systems  operate  on  the 
localizer  and  glide  path  receiver  output  signals  to  generate  roll  and  pitch 
attitude  commands  to  the  respective  autopilot  channels  to  force  the  aircraft 
to  fly  to  the  indicated  course/path  null.  Any  noise  in  the  indicated  track 
error  signals  causes  resulting  dispersions  in  the  actual  position  of  the  air¬ 
craft.  An  acceptable  automatic  ILS  approach  is  one  which  satisfies  pre¬ 
scribed  performance  criteria  such  as  control  column  activity,  aircraft  at¬ 
titude  and  displacement  dispersions  at  various  points  along  the  track  in  re¬ 
sponse  to  this  course/path  noise. 

The  present  figure  of  merit  for  localizer  course  and  glide  path 
acceptability  is  a  maximum  noise  amplitude  specification  which  does  not 
involve  approach  performance  criteria.  This  allowable  amplitude, 
which  is  reasonably  large,  is  an  outgrowth  of  many  past  compromises 
in  the  selection  of  practical  ILS  subsystem  component  tolerances.  These 
tolerances  are  adequate  for  the  present  use  of  the  ILS  to  the  200  foot  min¬ 
imum.  However,  if  lower  minima  and  eventually  final  landings  are  to  be 
accomplished  with  the  aid  of  the  ILS,  tighter  control  of  the  beam  charac¬ 
teristics  and  component  tolerances  should  be  provided.  Since  zero  noise 
amplitudes  cannot  practically  be  specified,  some  reasonable  value  of  a 
pertinent  measure  of  noise  must  be  used.  This  parameter  may  be  a  small¬ 
er  value  of  maximum  allowable  amplitude  based  on  standard  deviations 
which  provide  acceptable  tracking  performance;  or  it  may  be,  for  exam¬ 
ple,  a  parameter  which  additionally  specifies  the  course/path  noise  as  a 
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function  of  frequency.  The  latter  situation  was  assumed  in  Phase  I 
of  the  contract  work  statement,  which  indicates  that  the  bend  power 
spectral  density  is  a  probable  figure  of  merit.  The  Phase  I  analytical 
studies  were  designed  to  determine  the  extent  to  which  the  above 
parameters  may  be  used  as  figures  of  merit.  Section  Z  of  this  report 
contains  the  results  of  the  power  spectral  analyses  of  most  of  the  lo¬ 
calizer  and  glide  slope  data  taken  at  25  major  airports  in  the  U.  S. 
that  was  furnished  by  the  Federal  Aviation  Agency  under  Phase  I. 

This  FAA  data  was  derived  employing  theodolite- controlled  approaches 
wherein  the  aircraft  tracking  error  was  in  the  order  of  2.  5  feet,  the 
data  rate  was  approximately  1  reading  per  second;  and  the  ILS  course/ 
path  recorded  with  full  scale  recording  deflections  (+  2  inches)  of 
+  30  microamperes,  and  a  recording  speed  of  24  inches  per  minute. 

The  recorder  frequency  response  was  approximately  25  cps.  The  vis¬ 
ual  tracking  error  was  removed  from  the  test  recordings,  so  that  the 
resultant  tape  recordings  represented  the  actual  ILS  course  structure 
with  a  reasonable  degree  of  accuracy. 

The  effort  in  Phase  lA  included  the  following  studies; 

1,  It  determined  the  extent  of  the  usefulness  of  a  simple  filter 
response  technique  in  specifying  ILS  beam  acceptance.  This 
technique  involved  the  use  of  relatively  simple  filters  which  simu¬ 
lated  the  response  of  the  actual  airborne  systems.  The  trade-off 
between  filter  simplicity  and  correlation  with  the  full  scale  analog 
computer  responses  was  investigated. 

2.  It  determined  the  extent  of  the  usefulness  of  the  noise  vari¬ 
ance  program  technique  in  specifying  ILS  course/path  acceptance. 
This  technique  involved  the  use  of  suitably  shaped  random  noise 
signals  whose  total  power  or  variance  was  modified  as  a  function 
of  range  to  the  ILS  transmitters  to  simulate  the  effects  of  a 
course/path  noise  amplitude  specification.  Several  variance  pro¬ 
grams  were  studied  to  determine  the  maximum  allowable  course/ 
path  noise  which  could  be  tolerated  by  the  test  systems  for  low  ap¬ 
proach  operations. 
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The  results  of  the  stationarity  and  normality  tests  are  presented 
in  Section  2.  1  along  with  the  time  histories  of  the  course  noise  at  each 
facility.  The  power  spectra  are  presented  in  Section  2.  2.  Section  2.  3 
presents  the  closed  loop  simulated  time  responses  of  the  KC-135,  DC-7, 
and  TF-102  aircraft  under  the  actual  course  noise  conditions.  The  time 
responses  are  shown  for  the  two  lateral  control  systems  using  heading 
and  derived  beam  rate  as  track  damping  terms  in  the  approach  couplers. 
Section  2.  4  presents  the  results  of  attempts  at  correlating  the  actual 
track  responses  with  several  possible  course  power  noise  figures  of  mer¬ 
it.  The  conclusions  for  the  use  of  the  power  spectrum  as  a  figure  of  mer¬ 
it  for  approach  performance  are  given  in  Section  2.  4,  Section  2.  5  contains 
the  results  of  the  glide  path  stationarity  tests  for  some  of  the  facilities. 
Section  2.  6  contains  some  glide  slope  dispersion  data  for  the  KC-135, 

DC-7,  and  TF-102A  aircraft.  Section  3.  1  discusses  the  results  of  the 
simple  filter  technique  as  applied  to  the  three  test  aircraft.  Section  3.  2 
contains  the  results  of  the  course/path  noise  variance  program  technique. 
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Figure  1  Airport  Approach  Geometry 
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1.  PHASE  I  BEAM  ACCEPTANCE  STUDIES 

1.  1  LOCALIZER  STATIONARITY  AND  PROBMULITY  DISTRIBUTION  TESTS 


Figure  1  and  Table  I  indicate  the  geomi  trical  relationships  between 
the  localizer  and  glide  path  transmitters,  th..-  inarker  beacons,  and  the  run¬ 
way  at  each  of  the  test  facilities. 
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Figure  Z  shows  the  theodolite  corrected  course  noise  at  each 
of  the  localizer  facilities.  It  is  this  data  that  was  tape  recorded  and 
processed  to  provide  the  statistical  figures  of  merit.  As  indicated 
in  Reference  1,  knowledge  of  the  stationarity  and  normality  of  the 
course/path  noise  is  necessary  prior  to  deriving  the  power  spectra. 
As  an  example  of  the  procedures  outlined  in  that  report,  Table  II  in¬ 
dicates,  in  discrete  fashion,  the  time  history  of  the  variance  of  the 
Miami  Localizer.  Although  there  are  deviations  from  the  average 
variance,  they  are  not  considered  significant  and  the  course  noise  at 
this  facility  is  said  to  be  at  least  weakly  self- stationary.  Table  III 
indicates  the  results  of  the  other  localizer  facilities. 
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A  "good"  facility  is  one  for  which  none  of  the  subrecord  variances 
sii'Mificantly  exceeded  the  range  of  values  prescribed  iiv  tlie  Chi-square 
lest.  A  "fair"  facility  is  one  for  which  no  more  tli  tn  tvo  of  the  ten  sub¬ 
record  variances  significantly  exceed  the  pi  «'s«;  nln  d  values.  A  "poor" 
laeility  is  one  for  which  more  than  two  of  the  ten  subrecord  variances 
SI  ‘Ml ficantly  exceeded  the  prescribed  values.  The  results  indicate  that 
Hilly  the  St.  Louis  facility  can  be  considered  pi>t»rly  self- stationary.  The 
"Illy  significance  of  this  statement  is  that  a  single  power  spectrum  for 
'his  facility  would  not  be  as  useful  a  possible  figure  of  merit  as  it  would 
'  for  the  spectra  of  facilities  which  are  at  least  weakly  self- stationary. 

II  does  not  indicate  the  quality  of  the  course  noise  as  it  relates  to  air- 
<  raft  performance. 


TABLF:  III 

RESULTS  OF  STATION AHITY  TESTS 
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Figure  3  indicates  the  result  of  the  normality  tests  with  the  course 
noise  of  the  Miami  Localizer  facility;  and  Table  IV  is  a  tabulation  of  the 
mean  and  standard  deviations  for  this  and  fht'  other  facilities.  The  results 
of  these  tests  for  a  majority  ot  the  locali/.er  facilities  indicate  that  the 
course  noise  is  reasonably  Gaussian  in  character,  at  least  within  one  or 
two  standard  deviations.  These  results  imply  that  those  facilities  which 
were  weakly  self- stationary  are  probably  also  strongly  self- stationary  and 
that  a  single  power  spectrum  can  be  derived  in  most  i  ases. 
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l.l  LOCALIZER  POWER  SPECTRAL  TESTS 


Figure  4  indicates  the  results  of  the  power  spectral  tests  for  a 
nnajority  of  the  localizer  facilities.  A  general  trend  of  increasing  power 
with  decreasing  frequency  is  quite  evident.  Table  IV  also  includes  a  com¬ 
pilation  of  the  average  power  in  each  of  three  spectral  bands  as  well  as 
the  total  course  noise  power,  the  mean  power,  the  average  power  at  the 
middle  marker  and  the  average  power  between  the  outer  and  middle  mark¬ 
ers.  The  starred  facilities  had  directional  localizers  at  the  time  the  flight 
recordings  were  made.  The  NAFEC  and  Duluth  facilities  have  very  low 
noise  power  in  the  measured  spectral  bands.  The  total  random  noise  power 
at  the  NAFEC  facility  is  not  significantly  lower  than  the  other  facilities 
however,  and  this  indicates  that  considerable  power  is  located  in  the  very 
low  frequency  band  of  the  power  spectrum  which  could  not  be  measured  di¬ 
rectly.  A  glance  at  the  time  history  for  this  facility,  in  Figure  2,  indicates 
that  this  is  true.  Unh-rtunately,  as  noted  in  Reference  7,  power  measure¬ 
ments  at  these  very  low  frequencies  are  unreliable  due  to  the  tape  splicing 
effect  and  bias  errors.  As  noted  in  Table  IV  the  mean  course  noise  power 
level  at  some  localizer  facilities  is  quite  high.  It  cannot  be  determined 
whether  this  mean  deviation  is  due  to  a  theodolite  offset  or  is  actually  due 
to  a  misalignment  of  the  beam  structure.  The  effect  of  the  mean  course  de¬ 
viations  was  subtracted  from  the  results  of  the  aircraft  response  tests  since 
any  misalignments  can  be  accounted  for  in  the  ILS  commissioning  process, 
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Figure  4  Power  Spectrum  -  Localizer  Facility 
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2.  3  LOCALIZER  TIME  RESPONSE  TESTS 

Figures  5  and  6  are  block  diagrams  of  the  lateral  approach  coupler 
systems  that  were  simulated  on  the  analog  computer.  The  significant  dif¬ 
ference  between  the  two  systems  is  the  use  of  a  relative  heading  signal  as 
a  track  damping  term  in  Figure  3  in  lieu  of  the  derived  beam  rate  signal  in 
Figure  6. 


Figure  5  Heading  Damped  Coupler 


Figure  6  Beam  Damped  Coupler 


The  heading  damped  system  requires  a  beam  error  integration  sig¬ 
nal  to  provide  tight  control  under  crosswind  conditions.  The  rate  filter  con¬ 
stants  and  were  chosen  on  the  basis  of  track  stability  and  are  con¬ 
sistent  with  the  lowest  first  order  lag  break  frequency  which  provides  rea¬ 
sonable  course  noise  rejection.  Figures  7  and  8  show  the  simulated  time 
responses  of  the  KC-135  autopilot/air  craft  with  each  coupler  system  with 
no  course  noise.  The  rate  damped  system  track  response  is  somewhat 
faster  than  the  heading  damped  system  due  to  the  airplane's  natural  weath- 
ercocking  capability.  This  would  lead  one  to  suspect  that  the  track  devia¬ 
tions  under  noise  conditions  might  be  smaller  with  the  rate  damped  system 
at  the  expense  of  increased  roll  and  stick  activity  due  to  the  derived  rate 
signal.  This  point  is  illustrated  by  comparing  the  response  of  each  KC-135 
coupler  at  the  NAFEC  facility  in  Figures  9  and  10.  This  facility  has  a  117 
foot  waveguide  localizer  so  that  the  control  activity  with  either  coupler,  is 
relatively  low.  However,  the  rate  system  generates  two  to  three  times  more 
activity  than  the  heading  system.  It  is  also  worthwhile  to  note  the  difference 
between  the  actual  course  noise  and  the  indicated  course  error  in  Figure  9. 
The  aircraft  track  deviation  (y)  from  the  ideal  center  due  to  course  noise  is 
not  much  less  than  the  deviations  in  Figure  10  with  the  heading  system.  The 
result  is  that  although  the  rate  system  offers  improved  tracking  responses 
on  beams,  the  responses  are  significantly  degraded  when  the  ILS  signal  con¬ 
tains  course  noise. 
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Figure  7 

KC-135A  Track  Response 
Heading  Damped  Coupler 
No  Noise 


Figure  8 

KC-135A  Track  Response 
Rate  Damped  Coupler 
No  Noise 


Figure  9 

KC-135A  Track  Response 
Heading  Coupler 
NAFEC  Localizer 


Figure  10 

KC-135A  Track  Response 
Rate  Coupler 
NAFEC  Localizer 


Tables  V  and  VI  indicate  the  aircraft  responses  to  the  noise  at 
each  facility  at  the  middle  marker,  the  runway,  and  the  touchdown  re¬ 
gion  for  each  type  of  approach  coupler.  Figures  11  and  12  are  time 
histories  at  the  Miami  Localizer  facility  with  the  KC-135A  aircraft 
with  each  coupler  system. 


TABLE  V 

DISPERSION  DATA 
HEADING  SYSTEM  KC-135 


raiility 

Middle  Marker  Region 

y  i  ♦ 

Ft  Deg  Deg  Deg 

Deg 

Runway  Threshold 

y  4  • 

Ft  Deg  Deg  Deg 

‘’I. 

Deg 

mm 

rwy 

S.“ 

Fort  Worth 

<1 

U 

0.  2 

0 

0 

10 

tO.  1 

0 

0 

0 

5 

10 

.  2 

Detroit 

a 

tO.  S 

0 

0.  ) 

0.6 

10 

0 

0 

0 

0 

10 

10 

0 

Baltimore 

9 

0 

0 

0.  5 

0.  2 

0 

0 

0.  7 

0.  2 

0.  1 

5 

0 

0 

San  Francisco 

fO.  4 

0.  7 

0.  5 

1 

20 

♦  0.  2 

1 

0.  5 

1.4 

20 

15 

NAFCC 

< 

♦  0  5 

0.  75 

0.  5 

0. 

b 

0 

0.  5 

0.  1 

0.  1 

6 

8 

t 

Duluth 

H 

0 

0 

0 

0 

10 

♦  0.  1 

0 

0 

0 

8 

5 

.  2 

Miami 

57 

-0.  5 

0 

0.  75 

0.  5 

14 

♦  0.  5 

0.  5 

0.  75 

0.  5 

to 

20 

Chicago 

25 

»0  7 

0.  5 

0.7 

0.2 

15 

0 

0.2 

0 

0.  1 

20 

18 

.  45 

Uos  Angeles 

|5 

eO.  5 

i  H 

0.  42 

0.45 

♦  0.2 

0.  75 

0.  7 

0.  8 

12 

10 

.  9 

Dallas 

fc5 

eO.  5 

2 

0.9 

0.7 

20 

el 

0.8 

0.  5 

0.9 

to 

10 

. 

Idiewtld 

h 

eO.  1 

0.  2 

0 

0 

•  0  1 

0.  t 

0 

0 

4.  t 

4.  t 

0 

Minneapolis 

O 

fO.  5 

0 

0.  7 

0.5 

♦  0.25 

0.  5 

0.  ) 

0.  2 

to 

5 

.  5 

Houston 

10 

0 

0.  8 

0.® 

0.  2 

15 

♦  0.  5 

0 

0 

0.  1 

15 

Washington.  O  C. 

5 

♦  0.  1 

1 

0  0 

0  7 

0 

♦  0.  4 

1.2 

0.  7 

1 

D 

1 

. 

St.  Louis 

to 

eO  5 

0.  H 

0.  7 

0.  5 

15 

0 

0 

0.  7 

1 

15 

10 

0 

Cleveland 

40 

*0.  1 

0.2 

0.  45 

0.  5 

20 

♦  0.  8 

0.2 

0 

0.  t 

to 

15 

. 

Ontiriv* 

20 

0.  1 

2 

1.2> 

0.  i 

to 

♦  0.  > 

0.  •> 

0.  4 

0.  » 

25 

to 

Louisville 

b 

rO.  t> 

0.  4 

0.  5 

0.  5 

♦  0  5 

0.  t 

0 

0.  5 

7 

New  Orleans 

DO 

»0.7 

0.  5 

0.  t 

0.  t 

»2 

♦  0.  » 

0.2 

0 

0.  4 

45 

28 

Bi  rmingham 

45 

•  0  t 

0.  5 

0.4 

0.4 

45 

♦  0.  7 

0 

0.  5 

0.  7 

49 

to 

Burbank 

90 

•  2.  5 

4 

1.4 

0  115 

12 

•  2.0 

4 

t 

1.4 

tt- 

Atlanta 

12 

•  0-  5 

t.  5 

2.  h 

I.  2 

♦  1.5 

1 

0.  5 

1 

20 

Note  •  <.lu«ing  to  runway  <.enterlino 

•  V  •  departing  runway  crntt-rlin«*  j  ) 


Figure  11  KC-135A  -  Heading  Coupler  -  Miami  Localizer 
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This  facility  was  chosen  as  being  representative  of  nondirectional 
localizers.  The  advantage  of  the  NAFEC  installation  is  apparent  when  com¬ 
paring  Figures  9  and  10  with  11  and  12  since  the  attitude  and  track  deviations 
are  significantly  improved  with  the  lower  noise  content  of  the  directional  lo¬ 
calizer. 


TABLE  VI 
DISPERSION  DATA 
RATE  SYSTEM  KC-135 


Facility 

Middle  Marker  Region 

Y  t  <76^ 

ft  deg  deg  deg  deg 

Runway  Threshold  Region 

Y  *  ♦  0^  06, 

ft  deg  deg  deg  deg 

«y 

mm 

rwy 

ft 

Touch¬ 

down 

Region 

(y  ck 

rwy 

TD  TD 
ft  deg 

Fort  Worth 

IS 

-0.6 

2.  S 

1.4 

1 

19 

0 

0.  8 

0.  8 

1.3 

12 

10 

0 

Dotroit 

IS 

+1.S 

1 

2.4 

2.7 

30 

•  0.  8 

1 

1.4 

1.6 

IS 

18 

1.  8 

Baltimore 

0 

•1-2 

0 

1.4 

2 

20 

♦  1 

2 

2.  4 

1.6 

18 

6 

.9 

San  Francisco 

6 

•1.3 

1.2 

2.  2 

4.9 

21 

+  1.8 

2.  8 

2.  7 

7.  1 

18 

10 

. 

NAFEC 

20 

0.6 

1.  s 

1.  7 

0.  S 

0 

0 

0.  4 

0.  6 

0.  3 

13 

13 

.  8 

Duluth 

IS 

0.2 

0 

0 

0.  7 

6 

0 

0.  3 

0 

0.  7 

7 

9 

1.3 

Miami 

23 

-2.  3 

2 

2.2 

3 

7 

-1 

1 

2.  1 

3 

10 

10 

0.  8 

Chicago 

27 

-1.  s 

4 

2.4 

1.9 

21 

0 

1 

1.4 

1.3 

22 

22 

3.  8 

Loa  Angeles 

4 

0 

3.  2 

2.  8 

2.6 

6 

0 

1.2 

2.  S 

2.8 

6 

6 

2.6 

Dallas 

3S 

to.  3 

6.  S 

4.4 

3.  1 

10 

+  0.  1 

0 

4 

S 

28 

8 

1.2 

Idlewild 

3 

to.  3 

2 

2 

1.  3 

3 

-0.2 

0 

0.  8 

1.2 

3 

4 

.4 

Minneapolis 

6 

tO.  4 

0.  6 

0.9 

1.3 

K 

+  0.  6 

1.6 

1.7 

0.  8 

10 

18 

1.4 

Houston 

0 

0 

0 

1.  4 

0.  s 

S 

0 

.  S 

.  8 

0 

8 

. 

. 

Washington.  0.  C 

4S 

♦1.  3 

1.8 

0.  8 

2.6 

18 

+0.2 

2 

3.  2 

4.  8 

20 

18 

3.7 

St.  Louis 

7.8 

t0.2S  2 

8 

3 

18 

-1.2 

1 

2 

7 

18 

20 

. 

Cleveland 

IS 

+1.4 

0.  8 

1.4 

2.9 

7 

♦  1 

0.  8 

1 

1.9 

7 

12 

2.8 

Ontario 

23 

♦1.3 

3.  8 

4 

1.3 

10 

-0.  3 

0.  8 

1.  3 

2.  4 

16 

11 

1.8 

Louisville 

30 

-1.78  6 

3.  8 

4.  8 

24 

♦  1.  8 

0 

1.8 

3.  1 

40 

. 

, 

New  Orleans 

20 

+1.2 

2 

1.2 

2. 1 

30 

.  8 

3 

1.4 

2.6 

18 

, 

1.  3 

Birmingham 

33 

-0.4 

2 

1.7 

2.2 

3 

+  1.8 

2 

1.4 

3.  3 

33 

8 

1.  1 

Burbank 

120 

♦1.3 

*Um 

8 

1 

24 

-0.7 

8 

4.8 

2 

38 

. 

, 

Atlanta 

110 

+1 

8 

4.4 

1.8 

3 

+2.  4 

2 

1.4 

4 

70 

• 

• 

Figure  12  KC-135A  Track  Response  -  Rate  Coupler 

Miami  Localizer 
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The  localizer  course  noise  at  some  of  the  facilities  was  also 
applied  with  each  coupler  system  to  the  TF-102A  and  DC-7C  aircraft. 
Table  VII  contains  the  dispersion  data  for  the  TF-102  and  Table  VIII 
contains  the  data  for  the  OC-7C.  Comparing  these  results  with  those 
shown  in  Tables  V  and  VI  for  the  KC-135A  indicates  that  they  are  in 
reasonable  agreement. 


TABLE  VII 

TF-102A  AIRCRAFT  LOCALIZER 
DISPERSIONS  AT  MIDDLE  MARKER  AND  RUNWAY 


HEADING  DAMPED  SYSTEM 


Facility 

Y 

Middle  Marker 

*^deg  ^deg  ^♦deg 

Y 

Runway  Threshold 

'^deg  ^deg  *^^deg 

mm 

-rwy 

ft 

Touch¬ 

down 

-TD  TD 
ft  del 

Rmarla 

Nalec 

T~ 

t,.25 

1 

rs 

5 

1 

75 

75 

8 

St.  Louie 

22 

+  .  4 

.  8 

11 

0 

1 

.9 

15 

10 

.  5 

Burbank* 

50 

-.4 

5 

8 

+.  1 

- 

- 

30 

- 

- 

San  FrancUco 

2 

t.  2 

1.2 

8 

- 

0 

1.6 

s 

• 

- 

Dallas 

40 

+  .  3 

.9 

5 

• 

1.2 

1.2 

25 

10 

1.  I 

BEAM  RATE  DAMPED  SYSTEM 

Naiac 

15 

+  1.2 

3.8 

3.  5 

0 

2 

2.  6 

2.  5 

10 

5 

2.4 

St.  Louis 

0 

+  1 

3 

5.  8 

12 

tl 

0 

3 

10 

18 

5.  5 

Burbank* 

150 

+  .75 

10 

10 

24 

- 

- 

- 

45 

• 

- 

roll  lint. 

Dallas 

15 

-1 

7 

7 

16 

-1 

3 

5 

15 

10 

3.  7 

iting  iiv 

sidenm 

*See  Tabic  VIU 

not« 

_ 

-  -  — .  — ^ 

TABLE  VIII 

DC -7  AIRCRAFT  LOCALIZER 
DISPERSIONS  AT  MIDDLE  MARKER  AND  RUNWAY 


HIAMHO  BAUPSa  CgVPfcim 


Facility 

Y 

Mldeie  Marker 

■’tdei 

Y 

RSnway  Threahold 

*deg  ^tdeg  ®*deg 

mm 

rwy 

ft 

pooch, 
down  ^ 

«r  "♦ 

rwy  rwy 
TO  TO 
ft  deg 

Nafcc 

2S 

..25 

.  5 

.  2 

10 

4.25 

.2 

.  3 

15 

8 

.45 

St.  Louis 

21 

♦  .8 

.  45 

.  82 

15 

..  1 

.  42 

1 

5 

18 

.45 

Burbank* 

M 

■1.6 

2.2 

2 

12 

♦  .  6 

.  52 

.  42 

40 

. 

. 

San  Francisco 

6 

f.  2 

.62 

.86 

15 

0 

.  7 

1.5 

10 

. 

. 

Osllaa 

40 

♦  .  5 

.91 

1 

2 

♦.8 

.  55 

1.  1 

30 

4 

.6 

BEAM  RATE 

DAMPED  COUPLER 

Nafac 

30 

f  1 

2.2 

1 

15 

-.4 

1.6 

1.4 

15 

10 

.7 

St.  Louis 

30 

-.4 

2.  5 

3.  5 

22 

..  6 

1.0 

5 

27 

20 

Burbank* 

120 

.3.  4 

4 

4 

12 

•  1 

3.  1 

4 

50 

. 

. 

San  Francisco 

24 

♦  .  7 

1.  4 

4.  4 

10 

♦  1 

2.2 

7 

10 

. 

. 

Etallas 

14 

0 

3.0 

2 

25 

0 

1.  5 

1.  4 

10 

15 

1 

'Burbank  local 

1 

Bcr  tranamittcr  is  900  feet  in  1 

ront  of  the  approach  end  of  the  ILS  ru 

nwey 

1  Tha  runway  parameters  for  this  facility  were  actually  taken  700  feet  in  front  of  the  trana* 

1  mitter. 
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The  conclusions  in  Section  2.  4  therefore  also  pertain  to  these 
aircraft. 

The  95%  probability  levels  of  course  noise  in  the  middle  marker 
runway  region  are  listed  in  Table  IV  for  each  facility.  They  indicate  that 
the  test  facilities  have  95%  noise  amplitudes  between  3  and  20  ^a.  How¬ 
ever  only  nine  of  the  selected  facilities  would  pass  a  5  p.a  95%  noise  am¬ 
plitude  specification  while  19  facilities  would  pass  a  10  pa  specification. 
It  should  be  noted  that  all  the  facilities  pass  the  present  15  pa  specifica¬ 
tion  in  this  region  of  the  approach,  except  Burbank. 

The  airborne  system  configurations  for  these  aircraft  are  identi¬ 
cal  to  those  used  with  the  KC-135A  aircraft.  Different  gain  parameters 
are  used  however  to  provide  adequate  body  axis  and  track  stability. 

TABLE  IX 

LATERAL  CONTROL  SYSTEM  PARAMETERS 


Autopilot 

oil  and  rud  aarvo  natural  fraq  (rad/aae) 

KC>)3SA 

10 

DC-7C 

10 

Tr-102A 

10 

ail  and  rud  aarvo  damping 

.7 

.  7 

flrat  ordar  lag 

roll  att  rata  gain  -  (dag  ailaron/dag  roll  rata) 

l.l 

1.1 

9 

roll  att  gain  •  (dag  ailaron/dag  roll) 

2.  1 

2.1 

1.79 

yaw  rata  gain  •  (dag  ruddar/dag  yaw  rata) 

9.S 

9.9 

2.9 

coordination  gain  •  (dag  ruddar/dag  roll) 

1.99 

1.9 

.29 

roil  att  limit  •  ig  (dog  roll) 

10 

10 

6 

Haading  Couplar 

tiiaplacamant- gain  Ko  (dag  roll/dog  boam  orror) 

14.9 

14.9 

29 

Intogral  gain  -  K|  (dag  roll/aae^lag  baam  arror) 

.  149 

.  149 

.40 

haadtag  gain  •  (dog  roll/dog  hdg  orror) 

.99 

.  99 

2 

doaonaitiaatton  ratio  •  R 

9:  1 

9;  1 

9;  1 

Rato  Couplar 

Jlaplacamant  gain  •  Kq  (dog  roll/dog  boam  arror) 

7.9 

19 

19.9 

•  Kr  (dag  rol^ag/aac  baam  orror) 

190 

190 

910 

rota  tima  conataat  •  r|^  (aoeonda) 

2.9 

2.9 

2.2 

daaanaltlaation  ratio 

9:  1 

9:  I. 

9;  1 

raeaWor  lag  •  r  ( aaconda) 

.  9 

.  9 

.  9 
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Table  IX  provides  a  listing  of  the  pertinent  automatic  control 
system  parameters  that  were  used  with  each  aircraft-coupler  combina¬ 
tion. 

Table  X  presents  a  compilation  of  the  lateral  aerodynamic  para¬ 
meters  used  in  the  study. 

TABLE  X 

LATERAL  AERODYNAMIC  PARAMETERS 


Atrtpsed  (ft/tec) 

KC-135A 

210 

DC-7C 

202 

TF-102A 

290 

*»c 

2.611  X  10^ 

1. 32  X  10* 

.0192  X 10* 

4.59  X  10* 

1.  57  X  10* 

.  144  X 10* 

•lug -11^ 

0 

0 

-.0254x  10' 

W 

pounds 

160,000 

100,000 

23,400 

r*d*‘ 

..004 

-.086 

-854 

r«d*‘ 

*  .  00025 

*  .0163 

*  . 0224 

'-P 

rad** 

rad** 

-.385 

-.  531 

-.216 

% 

♦  .  33 

4  .  358 

*  .252 

rad** 

-.215 

-.  0774 

-.120 

C_ 

rad** 

-.0017 

-.  1604 

-.  0401 

"6R 

rad  ^ 

rad'* 

-.00016 

-.0152 

t  . 0344 

-.17 

-  . 1604 

-.291 

S 

rad‘* 

-.18 

-.118 

-.0912 

rad** 

*  .1346 

f  . 1 346 

4  . 12 

rad** 

-.687 

-.  676 

-.  561 

rad** 

*  .41 

0 

*■  .831 

r 
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L.  4  CORRELATION  OF  LOCALIZER  TIME  RESPONSES  AND  NOISE 
POWER  DATA 

Figures  13  through  21  are  attempts  at  correlating  the  output  time 
responses  in  Tables  V  and  VI  at  the  various  facilities  with  the  respective 
course/path  noise  power  parameters  in  Table  IV.  Figures  13,  14,  and 
15  are  graphs  of  lateral  deviation,  RMS  roll  attitude  and  RMS  aileron  de> 
flection  at  the  middle  marker  and  the  runway  of  a  majority  of  the  localizer 
facilities  with  respect  to  the  average  power  in  the  course  noise  spectral 
bands. 


A  f/At7  Atr  /TiyA/H/^r  (  A'A9rr  .fYSTfA4 


Figure  13  Lateral  Dispersions  vs  Average  Noise  Power  in  Spectral  Bands 
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The  low  frequency  band  Is  defined  from  .  02  to  .  1  cps,  the  mid¬ 
frequency  band  from  .1  to  .6  cps,  and  the  high  frequency  band  from  .  6 
to  2  cps.  The  rate  damped  system  does  provide  somewhat  smaller  lat¬ 
eral  dispersions  than  the  heading  system  but  the  roll  and  aileron  activity 
is  nearly  four  times  the  same  dispersions  for  the  heading  system. 

- -  o^peo  sysre/if  /tr  At/ooLt 

- aysnsM  yfr  Af/ooLe  KlflAHSA 

a  a  ACAL  /trr/ruoe  Ar  AHJMvmy- He/mifA}  i  ffme  sysreM 


Figure  14  RMS  Roll  Attitude  Dispersions  vs  Average  Noise  Power  in  the 

Spectral  Bands 
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There  are  weak  trends  toward  increasing  aileron  and  roll  activity 
with  average  spectral  course  noise  power  in  the  mid  and  high  frequencies. 
However,  there  does  not  appear  to  be  any  significant  trend  of  this  kind  for 
lateral  dispersion  at  any  frequency.  This  is  due  in  part  to  the  inherent 
problem  of  attempting  to  correlate  an  instantaneous  output  parameter  with 
an  average  lumped  course  noise  parameter  such  as  power. 


- tfOMreio ftr  ivfiuoLe 

pmrt  M'srrM  nr  tt/oote 

A  A  ntLEttoN  po-itr /ON  nr  i^nn-nr-ngnoiN^  t pnre  SYsreM 


Figure  15  RMS  Aileron  Position  Dispersion  vs  Average  Noise  Power  in 

Spectral  Bands 
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Figure  16  shows  the  relationship  of  the  three  output  functions  with 
the  total  random  course  noise  power  as  denoted  by  the  average  variance 
of  the  course  noise  measured  during  the  normality  tests  for  each  facility. 
Here  again  there  are  only  weak  trends  toward  increasing  roll  and  surface 
artivitv  with  increasing  total  course  noise  power. 
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Figure  17  shows  the  relationships  of  the  output  parameters  to 
the  mean  value  of  the  localizer  course  noise.  As  expected,  only  the 
lateral  deviation  is  reasonably  correlated  with  this  noise  parameter. 
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Figure  17  Aircraft  Dispersions  vs  Mean  Localizer  Beam  Error 
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^  18  shows  the  output  parameters  as  functions  of  mean  noise 
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Figure  18  ,A  , .  -  raft  Dispersions  v  ,  Total  Noise  Power 


Figures  19  and  20,  respectively,  show  the  variation  of  the  output 
parameters  with  the  average  random  course  noise  power  in  the  region  of 
the  middle  marker  and  the  region  between  the  middle  and  outer  markers. 
The  weak  correlations  of  the  output  parameters  in  Figure  15  are  not  ma¬ 
terially  enhanced. 
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Figure  19  Aircraft  Dispersion  vs  Average  Random  Power  Near  Middle  Marker 
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Figure  21  shows  the  relationship  between  the  total  random 
power  and  the  mean  power  and  the  average  spectral  power  in  the 
three  bands.  The  lack  of  correlation  is  quite  apparent. 
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Figure  21  Average  Spectral  Noise  Power  vs  Total  Noise  Power 


2.  5  GLIDE  PATH  STATIONARITY  TESTS 


Figure  22  shows  the  actual  theodolite  corrected  glide  path  noise 
at  each  of  the  test  facilities.  Refer  to  Figure  1  and  Table  I  for  the  geo¬ 
metrical  relationships  between  the  glide  slope  transmitter,  the  markers, 
and  the  runway. 

The  glide  path  noise  at  some  of  the  facilities  was  tested  for  sta- 
tionarity.  All  of  the  selected  sites  failed  this  test  due  primarily  to  sig¬ 
nificant  shifts  in  the  mean  and  standard  deviations  of  the  noise  subrecords 
near  the  runway.  Table  XI  shows  the  variance  fluctuations  for  five  fa¬ 
cilities.  It  is  impossible  to  determine  from  the  records  whether  this  path 
noise  characteristic  in  the  latter  portion  of  the  approach  actuall  exists  or 
is  apparent  path  noise  due  to  aircraft  motion,  low  theodolite  data  rate,  or 
inaccuracies  in  the  manual  theodolite  technique  at  close  range.  In  any 
event  the  records  cannot  be  used  to  derive  power  spectral  data.  However, 
the  alternate  techniques  discussed  in  Section  3  are  applicable  to  the  glide 
path  problem. 
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TABLE  XII 

LONGITUDINAL  AERODYNAMIC  PARAMETERS 
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TABLE  XIII 

LONGITUDINAL  CONTROL  SYSTEM  PARAMETERS 
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2.  6  GLIDE  PATH  TIME  RESPONSE  TESTS 


Figure  23  Indicates  the  control  systctn  that  was  used  to  evaluate 
the  effects  of  the  glide  path  noise.  The  approach  coupler  is  of  a  standard 
form  using  path  error  displacement  and  integral  terms  to  generate  a  pitch 
command.  The  autopilot  pitch  attitude  control  system  is  used  to  damp  the 
glide  slope  track  response.  Tables  XII  and  XIll  are  compilations  of  the 
aerodynamic  and  control  system  parameters,  respectively,  that  were  used 
with  each  test  aircraft. 


Figure  23  Glide  Path  Control  System 

The  simulated  aircraft  were  trimmed  to  fly  the  published  glide 
slope  angles  in  the  tests  which  involved  the  path  noise  at  each  facility. 
The  noise  was  injected  into  the  system  at  the  receiver  output  since  the 
original  flight  test  recordings  were  derived  from  that  source.  Para¬ 
meters  such  as  pitch  attitude,  range,  altitude  error  with  respect  to  the 
ideal  track  and  elevator  position  were  recorded  for  evaluation. 

The  autopilot  gain  parameters  were  adjusted  to  provide  adequate 
response  and  damping  to  step  signal  inputs.  Similar  tests  were  con¬ 
ducted  with  the  approach  coupler  parameters  to  provide  reasonable 
pulse  elevator  disturbances  at  several  close  range  positions  to  check  the 
effect  of  desensitization  on  system  stability. 
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TABLE  XIV 

GLIDE  PATH  RESPONSES  -  KC-135  AIRCRAFT 
ATTITUDE  DAMPED  TRACK  SYSTEM 
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TABLE  XV 

GLIDE  PATH  RESPONSES  -  DC-?  AIRCRAFT 
ATTITUDE  DAMPED  TRACK  SYSTEM 
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Tables  XIV,  XV,  and  XVI  are  compilations,  at  each  facility,  of 
the  performance  of  each  aircraft.  Figure  24,  as  a  typical  example,  is 
the  response  of  the  KC-135A  aircraft  to  the  Miami  Glide  Slope  signal. 


Figure  24  KC-135A  Track  Response  -  Miami  Glide  Slope 

These  results  indicate  that  the  pitch  dispersions  of  the  three 
aircraft  with  the  path  noise  do  not  differ  significantly.  The  autopilot 
configuration  tends  to  mask  the  differences  in  the  airframe  responses 
of  the  test  aircraft.  However,  the  altitude  deviations  at  the  runway 
range  from  two  to  ten  feet.  This  variation  is  significant  since  the  air¬ 
craft  is  only  fifty  feet  above  the  ground  at  the  runway  threshold.  If 
touchdown  dispersions  are  to  be  minimized,  the  glide  slope  should  pro¬ 
vide  more  precise  guidance  in  the  form  of  a  higher  signal  to  noise  ratio. 
The  alternate  path  acceptance  techniques  described  in  the  next  section 
were  used  to  determine  this  more  acceptable  level  of  path  noise. 

The  path  noise  figures  listed  in  column  one  of  Table  XIV  indi¬ 
cate  that  the  test  facilities  have  two- sigma  path  noise  values  which 
range  from  10  to  45  microamperes.  The  present  glide  slope  noise  spe¬ 
cification  does  not  define  an  acceptable  level  of  path  noise  in  the  region 
of  the  middle  marker  and  runway  although  30  ^a  is  taken  as  an  accept¬ 
able  value  prior  to  reaching  the  middle  marker.  Any  path  noise  ampli¬ 
tude  specification  should  at  least  reduce  this  amplitude  in  the  middle 
marker- runway  region.  The  95%  probability  level  of  noise,  in  column 
one  of  Table  XIV  indicates  that  only  two  of  the  selected  facilities  would 
fail  a  maximum  amplitude  noise  specification  of  30  pa;  and  nine  would 
fail  a  20  pa  specification. 


42-7709 
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3.  PHASE  lA  LOCALIZER  COURSE/GLIDE  PATH  ACCEPTANCE  STUDIES 


As  discussed  in  the  previous  section,  the  course/path  noise  power 
spectral  technique,  or  any  of  the  long  term  average  noise  power  estimates, 
have  not  allowed  the  approach  performance  to  be  specified.  This  lack  of 
statistical  correlation  is  due  in  part  to  the  finite  and  short  process  time  of 
the  approach.  These  results  indicate  however  that  there  is  some  improve¬ 
ment  in  correlating  the  average  output  dispersions  with  short  term  average 
course  noise  estimates.  Since  the  performance  of  the  landing  aircraft  is  of 
prime  importance  to  the  approach  process,  some  methods  of  including  the 
behavior  of  the  airborne  system  in  the  beam  acceptance  process  were  studied. 

The  two  techniques  that  were  studied  in  this  phase  of  the  program  were: 

1)  Course/Path  Noise  Variance  Program 

This  technique  provides  a  courso/^.'ath  noise  amplitude  specification 
as  an  acceptance  criterion  for  ILS  facilities.  It  operated  in  the  following 
way. 


A  random  noise  signal,  properly  filtered,  is  applied  to  the  simulated 
airframe/autopilot  coupler  combination  The  noise  variance  which  is 

a  measure  of  the  allowable  course  noise  excursion  is  arranged  to  vary 

in  a  smoothly  decreasing  manner  between  the  outer  marker  and  the  runway. 
For  example,  the  SC-79  maximvim  amplitude  recommendations  can  be  used. 

In  this  case,  the  localizer  variance  will  change  from  225  pa^  ‘ 

at  the  outer  marker,  to  6  pa^  *  5  pa)  at  the  middle  marker,  3500  feet 
from  the  runway. 

A  set  of  simulated  approaches  was  then  made  and  the  aircraft  track 
and  body  axis  dispersions  were  measured  at  the  runway  and  the  middle  mark¬ 
er.  If  the  2a  lateral  deviation  was  less  than  20  feet  and  the  2a  roll  attitude 
and  aileron  position  were  less  than  2  degrees,  the  variance  program  was  con¬ 
sidered  acceptable  The  actual  time  histories  of  the  responses  at  the  test  fa¬ 
cilities  using  the  theodolite  corrected  data  could  then  be  evaluated  to  deter¬ 
mine  the  individual  noise  variance  at  each  facility.  If  the  actual  variance 
changes  are  within  the  acceptable  range,  then  the  facility  can  be  used  for  low 
approaches.  The  converse  is  true  if  the  program  is  outside  the  acceptable 
range. 
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Several  ideal  programs  were  tested  to  determine  whether 
the  SC-79  amplitude  restrictions  are  realistic  or  are  too  severe. 

This  technique  has  the  advantage  of  simplicity  in  its  application  to 
the  present  FAA  flight  inspection  procedures.  However,  it  has  the 
disadvantage  of  not  explicitly  minimizing  the  false  acceptance  or 
rejection  probability  of  any  facility  for  low  approaches  since  the 
ideal  variance  prograon  is  artificially  produced. 

2)  Filter  Response  Techniques 

Technically  the  best  method,  short  of  actually  flying  a 
series  of  typical  aircraft  at  each  facility,  is  to  record  the  actual 
course/path  noise  and  present  it  to  a  complete  simulation  of  the 
airborne  system.  The  resulting  responses,  such  as  roll  and  pitch 
attitude,  surface  deflection,  and  track  deviation's  can  then  be  meas¬ 
ured  directly  to  determine  if  preselected  allowable  dispersions  have 
been  exceeded.  If  the  simulation  is  accurate  enough,  there  would  be 
very  little  uncertainty  between  the  analog  results  and  any  actual 
flight  test. 

However,  a  complete  analog  computer  facility  is  not 
readily  transportable  so  that  the  results  would  have  to  be  obtained 
at  some  central  FAA  computing  station.  A  set  of  simple  filters, 
which  could  be  carried  in  the  inspection  aircraft,  could  be  used  to 
generate  approximate  model  transfer  functions  of  the  desired  dis¬ 
persion  parameters.  The  noise  could  be  applied  to  these  circuits 
and  the  dispersions  determined  directly. 

3.  1  FILTER  TECHNICUE 

The  results  of  power  spectral  studies  in  Section  2. indicate 
that  it  is  the  individual  bends  of  the  course  noise  and  not  the  effebts 
of  any  long  term  average  of  the  noise  that  influence  the  performance 
of  the  approach.  Any  method  that  attempts  to  correlate  performance 
with  noise  should  account  for  these  short  term  course  noise  effects. 
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Figure  25  Pitch  Attitude  Longitudinal  Damped  Coupler  System 
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However  it  is  desirable  that  the  decision  to  commission  or  revali¬ 
date  an  ILS  facility  for  low  approach  use  should  be  made,  by  the 
flight  inspection  team,  on-site.  It  was  the  purpose  of  this  section 
to  determine  the  extent  to  which  the  rather  complicated  analog  com¬ 
puter  simulation  could  be  simplified  so  as  to  be  capable  of  being  air¬ 
borne  and  still  provide  useful  correlation  between  the  noise  and  air¬ 
craft  performance.  Simplifying  the  equivalent  filter,  that  represents 
the  ground  and  airborne  portions  of  the  ILS,  causes  some  uncertainty 
in  the  results  and  therefore  increases  the  probability  of  accepting  an 
unsatisfactory  beam  and  rejecting  a  useable  one.  It  was  felt,  however, 
that  any  simple  filter  should  not  exceed  third  or  at  most  fourth  order 
complexity  in  an  effort  to  maintain  simplicity  in  the  airborne  or  field 
environment.  The  analyses  were  conducted  using  the  KC-135  aircraft 
but  the  results  are  applicable  to  other  aircraft  in  the  approach  environ¬ 
ment.  In  the  lateral  channel,  the  filters  were  designed  to  provide  the 
roll  attitude  and  lateral  track  deviation  functions  with  respect  to  the 
localizer  course  noise  input.  In  the  pitch  channel  the  filters  were  de¬ 
signed  to  represent  the  pitch  attitude  and  vertical  track  deviations  with 
respect  to  the  glide  path  noise  input. 

3.  1.  1  Glide  Path  Filter  Response 

The  pitch  attitude  damped  coupler  system  shown  in  Figure  25 
was  used  as  a  model  for  the  glide  path  filters,  and  represents  the  com¬ 
plete  analog  computer  simulation  for  the  KC-135  aircraft.  It  is  based 
on  a  linear  two- degree -of-freedom,  constant  airspeed,  airframe,  and 
a  standard  pitch  attitude  referenced  autopilot.  The  system  gain  para¬ 
meters  are  listed  in  Tables  XII  and  XIII.  The  actual  transfer  functions 
of  pitch  attitude  to  beam  error  noise,  6/f(t),  and  altitude  to  beam  error 
noise,  h/f(t),  are: 


e  -5104  s  (s+.  025)(8-t-.  568)(B-t-2) _ 

K(x)f{t)  =  (g+.  03)(s+2.  02){82+1.  448+3.  31)(82+.  221  8+.  03)(82+17.  088+148.  6) 

h  _  -8790  (8+.025)(g+2) _ _ 

K(x)f(t)  (s+.03(s+2.02)(82+1.44  8+3.  31)(82+.  221  8+.  03)(82+l7.  08  8+148.6) 
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Figure  26  Bode  Plot  of  Pitch  Attitude/Beam  Responses 


TABLE  XVII 

SIMPLE  FILTER  PARAMETER  VALUES 


Lomttudtnal  Stmpt*  Ftlfr  L*ter»l  Simple  Filfr 


Parameter 

Third  Order 
Filter 

Fourth  Order 
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K 

T  Ti  C  «„ 
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2 

5 

.  8 

.  2 

2  2. 5  5  .8 

Roll  (6) 

112 

S  1 . 44  2.  24 

Altitude  (h) 

.06 

2 

S 
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The  factor,  K(x),  on  the  left  side  of  both  equations  represents 
the  desensitizing  factor  that  is  normally  used  in  the  approach  couplers 
to  account  for  beam  convergence.  A  frequency  response  of  the  function 
in  Equation  1  is  shown  in  Figure  26- 

The  simple  filter  design  was  a  third  order  approximation  whose 
parameters  were  adjusted  to  obtain  the  best  approximation,  for  a  step 
input,  between  the  exact  analog  computer  response  and  the  filter  response. 
The  filter  transfer  functions  used  were,  for  pitch  deviation. 


®F 

K{x)  f(  t) 


-KS 


(1+  TS)( 


~S2~ 

<»)_2 

F 


+ 


HE 

“F 


S  +  1) 


and  for  altitude  deviation. 


(3) 


hjr  -K 

- £ -  =  - ^ -  (4) 

K(x)  f(t)  (1+  TS)(-1-  X  g  ^ 

(l}p2  u,  p 

The  parameters  K,  T,  and  of  the  filters  were  optimized 
by  minimizing  the  difference 

I  ( 0  -  6p(K,  T,  ^  p,  ft;  p)  )=  minimum  (5) 

between  the  filter  (0p)  and  respective  computer  output  (6). 

Repeated  trials,  with  many  variations  of  the  filter  parameters, 
were  used  to  determine  the  best  approximation.  However  Equation  5  can 
never  be  equal  to  zero  since  the  computer  response  is  generated  by  an 
eighth  order  system  while  the  simple  filter  response  is  generated  by  the 
third  order  filter.  It  was  therefore  necessary  to  minimize  the  error  in 
Equation  5  in  a  particular  time  interval  of  the  step  response  tests.  The 
first  fifteen  seconds  was  taken  as  a  desirable  time  interval  since  it  tends 
to  match  most  of  the  short  term  attitude  and  track  motions.  The  results 
of  the  parameter  variations  are  shown  in  Table  XVU.  It  should  be  stressed 
that  the  optimization  was  based  on  responses  to  step  inputs.  Somewhat  dif¬ 
ferent  filter  parameters  would  perhaps  have  been  chosen  if  other  types  of 
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Figure  27  Ideal  Pitch  Attitude  and  Simple  Filter  Responses  to  Step  Noise 


6i-77H 
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input  signals  were  used  or  if  other  matching  time  intervals  were  desired. 
Figure  27  shows  the  analog  computer,  filter,  and  the  difference  signal 
representing  Equation  3  for  the  pitch  attitude  response.  The  initial  con¬ 
ditions  of  the  simulated  aircraft  in  this  test  series  were  such  that  it  was 
flying  along  a  2.  5  degree  glide  slope  and  the  results  in  Figure  27  show  the 
deviations  in  the  parameters  from  the  steady  state  flight  condition.  The 
primary  difference,  in  the  first  fifteen  seconds,  between  the  computer  and 
filter  output  signals  is  in  the  initial  responses  to  a  step  beam  displacement. 
Beyond  fifteen  seconds  the  error  signal  again  becomes  large.  These  differ¬ 
ences  are  due  to  the  compromise  in  high  and  low  frequency  matching  of  the 
simpler  filter  to  the  more  complex  computer  simulation.  In  order  to  reduce 
the  discrepancies  at  either  end  of  the  frequency  spectrum,  a  fourth  order 
filter  of  the  form 

Sf  -ks 

K(x)f(t)  "  (14-TS)(U  T|S)(_S1_  2^  s+i) 

tty 

was  also  implemented.  Figure  27  also  shows  the  response  of  this  filter 
after  an  identical  optimization  procedure  was  performed.  Table  XVIl  also 
lists  the  resulting  values  of  parameters  for  this  filter.  Comparing  the 
output  responses  of  each  filter  shows  that  no  significant  improvement  has 
been  obtained.  It  was  therefore  concluded  that  for  practical  considerations 
any  simple  filter  for  the  pitch  channel  should  not  exceed  third  order  since 
much  closer  correspondence  to  the  computer  output  can  only  be  accomplish¬ 
ed  when  the  filter  has  nearly  the  same  order,  and  therefore  complexity,  as 
the  computer  simulation. 

Figure  28  shows  the  response  of  the  third  order  altitude  filter  and 
the  computer  to  a  step  beam  displacement.  The  two  signals  are  matched 
to  almost  the  same  extent  as  the  pitch  responses.  The  altitude  filter  pa¬ 
rameters  are  also  given  in  Table  XVII. 

The  final  filter  designs  were  used  in  conjunction  with  the  complete 
simulation  to  compare  the  results  when  the  theodolite  corrected  course.Jioise 
was  applied  to  each. 

The  attitude  and  altitude  responses  at  the  runway  threshold  were 
tabulated  for  the  computer  simulation  and  the  filters.  If  the  filter  and  the 
computer  time  responses  coincided  throughout  the  tests  then  there  would 


43 


Figure  28  Ideal  Altitude  and  Simple  Filter  Response  to  Step  Noise 
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be  a  one-to-one  correspondence,  or  perfect  correlation.  A  plot  of  the 
two  output  signals  would  then  be  a  straight  line  with  a  45  degree  slope. 
Obviously  such  coincidence  in  the  time  domain  requires  that  there  be 
corresponding  coincidence  in  the  frequency  domain.  However,  compar¬ 
ing  the  actual  computer  transfer  functions  in  Equations  1  with  the  filter 
transfer  function  in  Equations  3  shows  that  only  some  degree  of  corre¬ 
lation  can  exist  between  the  two  signals.  Figure  26  includes  the  frequency 
response  of  the  simple  pitch  filter  for  comparison  with  the  computer 
frequency  response. 

A  qualitative  measure  of  the  degree  of  correspondence  between 
the  filter  and  the  computer  output  signals  is  the  correlation  coefficient, 
r.  The  correlation  coefficient  is  -l- 1  when  the  ideal  response  can  be 
exactly  predicted  from  the  simple  filter  response.  The  value  of  r  is 
zero  when  there  is  no  relationship  between  the  two  signals.  In  between 
0  and  f  1 ,  the  value  of  r  indicates  that  only  a  degree  of  predictability 
exists  and  that  some  error  must  always  be  expected.  The  sign  of  r  in¬ 
dicates  whether  the  simple  filter  response  increases  when  the  computer 
output  is  decreasing.  A  quantitative  measure  of  the  degree  of  corres¬ 
pondence  between  the  computer  and  filter  output  signals  is  provided  by 
the  correlation  coefficient  squared,  r^.  This  parameter  represents  the 
decimal  percentage  of  the  variance  of  the  simple  filter  response  that  can 
be  accounted  for  by  the  relationship  with  the  computer  response.  There¬ 
fore  if  the  correlation  coefficient  between  the  two  responses  is  0.  8  then 
64%  (r^  =  0.  64)  of  the  expected  variation  can  be  attributed  to  the  rela¬ 
tionship  between  the  responses.  The  remaining  36  per  cent  of  the  varia¬ 
tion  between  the  responses  is  due  to  other  factors  such  as  the  statistical 
nature  of  the  actual  beam  noise,  and  the  simple  filter  approximations  to 
the  complete  system  response. 

A  95%  variance  was  taken  as  an  acceptable  value  of  correlation 
in  an  effort  to  minimize  the  resulting  false  acceptance  probability  for 
both  the  glide  slope  and  localizer  filter  design.  Table  XVM  shows  a 
typical  calculation  for  the  correlation  coefficient  using  the  NAFEC 
Glide  Path  facility. 

Ten  glide  slope  facilities  were  used  to  test  the  degree  of  corres¬ 
pondence  between  the  filter  and  the  computer  responses.  Typical  com¬ 
puter  and  filter  responses  are  shown  in  Figure  29  using  the  NAFEC 
facility  as  an  example.  Table  XIX  is  a  tabulation  of  the  correlation 
coefficient  r  and  r^  for  the  responses  in  the  middle  marker /run way 
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Figure  29  Ideal  and  Simple  Third  Order  Filter  Responses  to  Beaun  Noise 

NAFEC  Glide  Slope  Facility 
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TABLE  XIX 

A  SUMMARY  OF  PITCH  ATTITUDE  CORRELATION  COEFFICIENTS 
FOR  TEN  GLIDE  SLOPE  FACILITIES 


Correlation 

Square  of  Correlation 

Glide  Slope  Facility 

Coefficient 

r 

Coefficient 
(Per  Cent) 
r^ 

LOUISVILLE 

-.  078 

0.  6 

BURBANK 

+  .  962 

92.  5 

HOUSTON 

+  .  953 

90.  8 

WASHINTON,  D.  C. 

+  .  673 

45.  3 

CHICAGO 

+  .  845 

71.4 

ONTARIO 

-.  228 

5.2 

ST.  LOUIS 

+  .  910 

82.  8 

SAN  FRANCISCO 

+  .  957 

91. 6 

NEW  ORLEANS 

+  .  833 

69.  4 

NAFEC 

+  .  982 

96.  4 

region  of  each  facility.  The  wide  range  of  values  of  r^  indicates  that 
there  is  significant  variation  in  the  correlation  between  the  filter  and 
the  computer  responses.  For  example,  the  NAFEC  facility  has  a  high 
degree  of  correlation,  (r^  =  .964)  while  the  Louisville  facility  has  a 
very  poor  relationship  (r^  =  .  006).  The  wide  range  of  r^  values  is  re¬ 
lated  to  the  uniformity  of  the  ILS  beam  error  signals.  As  shown  in  Sec¬ 
tion  Z  the  statistical  characteristics  of  the  glide  slope  signal  in  the  vic¬ 
inity  of  the  middle  marker  and  runway  are  not  uniform.  Therefore  the 
simple  filter  responses  fail  to  approximate  the  ideal  computer  responses 
with  equivalent  errors. 

The  frequency  responses  in  Figure  26  indicate  that  if  the  beam 
noise  has  significant  power  at  frequencies  below  .  05  cps  and  above  1  cps 
the  responses  will  be  poorly  correlated,  since  a  best  fit  occurs  only  in 
the  mid- frequency  band. 

Table  XX  summarizes  the  responses  of  the  filter  and  the  compute 
with  the  noise  at  each  facility  in  terms  of  their  means  and  standard  devia 
tions.  It  complements  the  lack  of  correlation  shown  in  Table  XIX. 
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3.  1.2  Localizer  Filter  Responses 

The  roll  attitude  and  localizer  track  deviation  filters  were  designed 
in  an  identical  manner  to  that  described  above  for  the  glide  slope  responses. 
The  system  that  was  used  is  shown  in  Figure  30.  The  simulator  parameters 
are  listed  in  Tables  IX  and  X. 


Figure  30  Heading  Damped  Localizer  Coupler  System 


Frequency  responses  of  each  lateral  parameter  are 
Figures  31  and  32. 


The  simple  filter  designs  were  third  order  systems 
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Figure  31  Bi>de  Plot  Roll  Attitude/Beam  Response 


Figure  33  Ideal  Roll  Attitude  and  Simple  Filter  Response  to  Step  Noise 
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Figure  34  Ideal  Lateral  Displacement  and  Simple  Filter  Response  to 

Step  Noise 
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The  optimization  of  the  filter  parameters,  K,  T,  provided  minimum 

error  between  the  computer  and  filter  signals  in  the  first  15  seconds  of  the 
response  to  step  localizer  beam  displacements.  Table  XVII  presents  the  re¬ 
sulting  parameters  and  Figures  33  and  34  show  the  responses  that  resulted 
from  the  step  beam  error  inputs.  As  in  the  case  of  the  glide  slope  responses, 
there  is  a  discrepancy  in  the  initial  transient  responses  and  during  the  long 
periods  after  15  seconds.  The  simulated  aircraft  was  initially  positioned  along 
the  runway  centerline  with  wings  level  and  the  displacements  due  to  the  ILS 
course  noise  are  deviations  about  these  initial  conditions. 

Ten  localizer  facilities  were  used  to  test  the  simple  filter  technique. 
Figures  35  and  36  show  the  typical  responses  at  the  New  Orleans  and  Chicago 
localizer  facilities.  The  roll  attitude  correlations,  r  and  r^  are  given  in 
Table  XXI.  The  resulting  variations  in  r^  indicate  that  the  degree  of  corre¬ 
lation  between  the  filter  and  computer  responses  is  not  adequate  (i.  e.  ,  less 
than  95%)  for  all  the  facilities.  Only  the  responses  using  the  noise  of  the  Bur¬ 
bank,  Chicago  and  Ontario  facilities  indicate  a  high  degree  of  correlation.  It 
is  interesting  to  note  that  there  is  improved  correlation  for  the  lateral  results 
over  the  longitudinal  results.  This  is  due  in  part  to  the  fact  that  the  localizer 
course  noise  as  shown  in  Section  2  is  more  statistically  uniform  or  stationary, 
than  the  glide  path  noise  in  the  vicinity  of  the  middle  marker  and  the  runway. 

TABLE  XXI 

A  SUMMARY  OF  ROLL  ATTITUDE  CORRELATION  COEFFICIENTS 
FOR  TEN  LOCALIZER  FACILITIES 


Localizer 

Correlation 

Square  of  Correlation 

Facility 

Coefficient 

Coefficient 

(r) 

(Per  Cent) 

(r2) 


LOUISVILLE 

.  944 

89.  1 

BURBANK 

.  998 

99.7 

BALTIMORE 

.  944 

89.  1 

WASHINGTON,  D.  C. 

.  907 

82.  3 

CHICAGO 

.995 

98.9 

ONTARIO 

.993 

98.7 

ST.  LOUIS 

.  642 

42.  1 

SAN  FRANCISCO 

.  813 

66.  1 

NEW  ORLEANS 

.928 

86.  1 

NAFEC 

.936 

87.6 
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A  summary  of  the  responses  of  the  exact  computer  simulation  and 
the  approximation  filters  is  listed  in  Table  XXII.  This  data  augments  the 
correlation  data  presented  in  Table  XXL 

The  results  of  the  simple  filter  technique  as  presented  above  indi¬ 
cate  that  there  is  not  sufficient  correlation  between  the  responses  of  a  rela¬ 
tively  simple  filter  and  a  complete  simulation  to  adequately  protect  against 
falsely  accepting  an  undesirable  facility  or  rejecting  a  desirable  facility.  To 
improve  the  degree  of  prediction  it  is  necessary  to  increase  the  complexity 
of  filter  design  to  almost  coincide  with  the  complete  simulation  and  to  ensure 
that  the  ILS  course/path  noise  characteristics  are  more  uniform  than  they 
are  at  present.  The  simple  filter  technique  is  therefore  not  recommended 
for  use  by  the  FAA  in  commissioning  ILS  facilities. 

TABLE  XXII 

A  SUMMARY  OF  LATERAL  STATISTICAL  CHARACTERISTICS 
FOR  TEN  LOCALIZER  FACILITIES 
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Figure  35  New  Orleans  Liocalizer  Facility  -  Filter  Kesponses 
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Figure  36  Chicago  L>ocalizer  Facility  -  Filter  Responses 
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3.  2  NOISE  VARIANCE  STUDY 


3.  2,  1  Localizer  Noise  Variance  Study 

The  object  of  this  portion  of  the  study  was  to  determine  whether 
an  amplitude  specification  can  be  used  as  an  acceptable  criterion  for  low 
approaches.  The  test  procedure  was  to  generate,  using  a  random  noise 
source,  a  signal  whose  standard  deviation  was  a  prescribed  function  of 
the  range  to  the  transmitter.  Figure  37  shows  the  assumed  noise  varia¬ 
tion  and  the  present  ILS  course/path  noise  amplitude  specification  assum¬ 
ing  that  the  maxumum  noise  amplitude  actually  represents  two  standard 
deviations.  The  maximum  course  noise  amplitude  at  the  outer  marker 
was  kept  at  a  2a  value  of  ^  30  pa  and  the  final  value  at  a  point  7000  feet 
from  the  runway  was  adjusted  during  the  tests  from  a  2o  value  of  5  pa  to 
30  pa.  The  course/path  noise  variation  between  the  outer  and  middle 
markers  was  chosen  to  be  a  uniformly  decreasing  function  instead  of  the 
present  step  specification  since  that  characteristic  is  more  likely  to  be 
encountered  in  practice.  The  course/path  noise  standard  deviation  was 
kept  constant  from  the  7000  foot  point  throughout  the  rest  of  the  approach. 


It 


Figure  37  Localizer  and  Glide  Path  Noise  Variance  Programs 

Since  the  simulated  course/path  noise  was  random  in  amplitude 
and  frequency  content  up  to  a  corner  frequency  of  3  radians/second,  ap¬ 
proximately  twenty  tests  were  conducted  for  each  condition  to  ascertain 
the  average  or  most  likely  responses. 
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In  the  lateral  study  the  pertinent  responses  were  considered  to 
be  the  lateral  deviation  from  the  ruim  iy  centerline  and  the  rolling  ac¬ 
tivity.  In  the  longitudinal  glide  slope  study,  the  pertinent  responses 
were  taken  to  be  the  altitude  deviation  from  the  ideal  glide  path  and  the 
pitch  attitude  activity. 

Table  XXIII  shows  the  results  of  the  lateral  study  and  Table 
XXIV  shows  the  longitudinal  results  for  the  KC-135,  TF-102,  and  DC-7 
aircraft. 


TABLE  XXIII 

SUMMARY  OF  LOCALIZER  RESULTS 
CONTINUOUS  VARIATION  OF  STANDARD  DEVIATION 
OF  SIMULATED  COURSE  NOISE 
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TABLE  XXIV 
SUMMARY  OF  RESULTS 

CONTINUOUS  VARIATION  OF  STANDARD  DEVIATION 
OF  SIMULATED  GLIDE  PATH  NOISE 
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In  an  effort  to  aid  in  validating  any  rourse/path  noise  amplitude 
specification  using  the  simulated  noise  source,  the  data  using  that  noise 
was  combined  with  the  results  of  similar  responses  using  the  actual  theo¬ 
dolite  corrected  ILS  recordings  as  described  in  Section  2. 
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Figure  38  Localizer  -  Lateral  Deviation  Results  -  Heading  System 


Figure  39  Localizer  -  Lateral  Deviation  Results  -  Beam  Rale  System 
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Figures  38  through  41  show  the  results  of  combining  the  actual 
and  simulated  course  noise  for  the  localizer  mode  of  operation.  The 
vertical  bars  are  capped  at  the  minimum  and  maximum  attitude  or  dis¬ 
placement  that  was  found  during  each  of  the  simulated  noise  test  sequences. 
The  single  marked  point  in  each  bar  represents  the  standard  deviation  of 
the  particular  parameter.  The  individual  points  shown  in  these  figures  rep¬ 
resent  the  results  of  the  response  tests  at  each  of  the  ILS  facilities  using  the 
theodolite  corrected  data. 
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Figure  40  Localizer  -  Roll  Attitude  Results  •  Heading  System 
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Figure  41  Localizer  -  Roll  Attitude  Results  -  Beam  Rate  System 
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The  results  of  the  simulated  versus  actual  course  noise  tests  only 
agree  with  one  another  on  the  average  since  the  simulated  noise  is  derived 
from  an  almost  perfect  random  noise  generator,  and  the  actual  course  bends 
are  not  nearly  as  statistically  stationary.  Another  reason  for  the  disagree¬ 
ment  is  the  fact  that  the  actual  course  noise  standard  deviations  do  not  change 
in  the  same  manner  as  the  simulated  noise  was  allowed  to  change.  However, 
for  both  types  of  localizer  coupler  systems,  a  two-sigma  noise  value  of  5  pa 
in  the  middle  marker- runway  region  would  restrict  the  expected  value  (mean 
plus  standard  deviation)  of  the  lateral  deviations  to  below  1  5  feet  and  the  ex¬ 
pected  value  of  roll  attitude  to  below  two  degrees.  It  is  believed  that  this  low 
value  of  course  noise  may  be  used  as  the  criterion  for  the  future  use  of  the 
ILS  localizer  since  the  additional  aircraft  responses  to  gusts  and  other  dis¬ 
turbances  would  increase  the  deviations.  Referring  to  Figures  38  through  41 
indicates  however  that  choosing  this  low  value  of  maximum  course  noise  forces 
the  rejection  of  some  of  the  present  localizer  systems  even  though,  in  the  case 
of  some  of  the  systems  with  higher  noise  content,  acceptable  performance  can 
in  fact  be  obtained. 

Even  more  importantly  there  is  one  facility  which  fails  the  performance 
specifications  but  passes  the  95%  probability  level  course  noise  amplitude  spe¬ 
cification.  Table  XXV  which  is  a  compilation  of  these  facilities  which  pass  or 
fail  the  performance  or  noise  specifications  shows  this  very  clearly.  A  con¬ 
flict  is  indicated  when  the  results  of  the  performance  or  noise  amplitude  spe¬ 
cifications  do  not  agree.  For  the  heading  damped  system  there  are  five  con¬ 
flicts  if  the  maximum  lateral  offset  specification  is  fifteen  feet.  The  results  at 
Miami,  however,  indicate  that  this  facility  would  pass  the  noise  specification 
and  not  pass  the  maximum  lateral  deviation  specification.  This  is  true  even  if 
that  specification  were  increased  to  20  feet  instead  of  15.  With  the  20  foot  spe¬ 
cification,  the  San  Francisco  and  Chicago  facilities  are  added  to  the  conflict 
list  wherein  the  noise  specification  is  failed  but  the  performance  is  acceptable. 

The  right  hand  position  of  Table  XXV  shows  the  same  type  of  results 
for  the  rate  damped  control  system.  In  this  case  the  conflicts  between  the 
performance  and  noise  specifications  also  include  the  undesirable  case  in  which 
the  Miami  facility'passes  the  amplitude  criteria;  but,  the  rolling  activity  is 
beyond  the  acceptable  2  degree  level.  All  the  other  conflicts  would  unneces¬ 
sarily  penalize  the  localizer  facilities.  Table  XXVI  shows  that  there  is  con¬ 
siderable  agreement  between  the  results  of  applying  the  performance  and 
noise  amplitude  specifications  to  the  three  test  aircraft.  Only  one  unsatisfac¬ 
tory  conflict  occurs  and  this  with  the  St.  Louis  facility,  with  the  DC-7  beam 
rate  coupler  system. 
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TABLE  XXV 

APPLICATION  OF  PERFORMANCE  AND  COURSE  NOISE  AMPLITUDE 
SPECIFICATIONS  TO  TEST  LOCALIZER  FACILITIES 
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NOTE: 

1 .  "Yes"  denotes  that  the  noise  and  its  effects  are  within  the 
specified  limits 

1.  "No"  denotes  that  the  noise  and  its  effects  are  not  within  the 
specified  limits 

3.  A  conflict  denotes  the  case  wherein  there  was  an  inconsistency 
between  the  noise  amplitude  and  performance  specifications. 
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TABLE  XXVI 

COMPARISON  OF  EFFECTS  OF  APPLYING  ACCEPTANCE  CRITERIA 
TOKC-135.  DC-7,  AND  TF- 102  AIRCRAFT 
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Yes  ^  Yes  Yes 

St.  Louis 

Yes 

Y  es 

No  (no) 

Yes 

Yes 

Burbank 

No 

No  (no) 

No  (no) 

No  (no) 

Dallas 

No 

No  (no) 

Yes 

Yes 

San  Francisco 

No 

Yes 

Y  es 

. 

Yes 

NOTE:  See  Note  for  Table  XXV 
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Figure  42  Glide  Path  -  Altitude  Deviation  Results  -  Middle  Marker 
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Figure  43  Glide  Path  -  Altitude  Deviation  Results  -  Runway 
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Figure  44  Glide  Path  -  Pitch  Attitude  Results  -  Middle  Marker 
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Figure  45  Glide  Path  -  Pitch  Attitude  Results  -  Runway 
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It  is  clear  that  this  approach  of  specifying  a  maximum  allowable 
noise  amplitude  is  not  the  best  technical  solution  to  the  problem  although 
it  does  provide  an  acceptable  solution  in  terms  of  immediate  and  relatively 
simple  implementation.  These  considerations  are  further  expounded  upon 
in  the  Appendix  where  the  problem  is  treated  using  statistical  decision  theory 
to  prove  that  a  complete  analog  simulation  offers  the  best  technical  approach 
to  specifying  ILS  beam  acceptance. 

3.  2.  2  Glide  Path  Noise  Variance  Study 

The  same  general  comments  that  were  just  made  for  the  localizer 
can  be  made  of  the  glide  path  results.  In  this  case,  as  shown  in  Figures  42 
through  45,  the  20  p.a,  95%  probability  path  noise  amplitude  in  the  middle 
marker-runway  region  is  acceptable  by  providing  pitch  attitude  standard  de¬ 
viations  of  less  than  one  degree  and  altitude  deviations  of  less  than  four  feet 
at  the  runway  threshold.  An  allowable  noise  amplitude  between  20  and  30  |xa 
could  be  chosen  for  those  coupler  systems  which  use  attitude  damping.  How¬ 
ever,  the  increased  disperions  of  those  coupler  systems  which  use  beam  rate 
damping  as  well  as  the  additional  responses  to  gusts  and  other  disturbances 
force  the  20  ^a  figure  to  be  recommended. 

Table  XXVII  is  a  compilation  of  the  glide  path  test  facilities  along 
with  information  as  to  which  passed  or  failed  the  various  performance  and 
path  noise  amplitude  specifications.  Four  conflicts  between  the  two  specifi¬ 
cations  are  noted  in  the  fourth  column  of  Table  XXVII.  Of  these  only  the 
Louisville  facility  has  the  undesirable  characteristic  of  passing  the  noise  am¬ 
plitude  specification  and  failing  the  performance  criteria.  Eight  of  the  test 
facilities,  including  Louisville,  pass  the  20  p.a  maximum  path  noise  amplitude 
specification.  Table  XXVUI  compares  the  results  of  applying  the  specifications 
to  the  three  test  aircraft  at  some  of  the  glide  path  facilities.  The  only  conflict 
occurs  at  the  St.  Louis  facility  which  would  be  rejected  unnecessarily  since 
the  performance  is  adequate  although  the  noise  amplitude  specifications  is  ex¬ 
ceeded. 


The  results  of  the  variance  programs  with  the  TF-102  system  were 
repeated  using  a  noise  bandwidth  of  1 . 96  rad/sec  in  addition  to  the  3.  14  rad/ 
sec  bandwidth  that  was  used  throughout  the  tests.  No  significant  differences 
were  found  to  exist,  as  shown  in  Figures  42  through  45. 

It  is  interesting  to  note  that  these  glide  slope  responses  were  derived 
uding  conventional  approach  coupler  systems  with  increased  desensitization 
rates  to  account  for  the  lower  altitude  approaches.  No  glide  slope  extension 
schemes  were  included.  It  would  seem  therefore  that  the  need  for  such  ex¬ 
tension  systems  would  decrease  if  the  present  ILS  facilities  could  be  guaran¬ 
teed  to  have  the  lower  noise  content. 
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TABI.E  XXVII 

APPLICATION  OF  PERFORMANCE  AND  PATH  NOISE  AMPLITUDE 
SPECIFICATIONS  TO  TEST  GLIDE  PATH  FACILITIES  KC-135 


Facility 

95%  Prob 
Noiae  Ampl 
^  20  ^a 
MM-Rwy 

^‘‘Rwy  ^  ^ 

®‘’Rwy  ^  2  Deg 

Conflict 

Detroit 

• 

San  Franciaco 

Yea 

Yea 

Atlantic  City 

Yea 

Yea 

Yea 

Miami 

Yea 

Yea 

Chicago 

Yea 

Yea 

Idlewild 

Yea 

Yea 

Yea 

St.  Lioula 

• 

• 

Ontario 

No 

Yes 

Minneapolia 

Yea 

Yes 

Houaton 

Yea 

Yea 

Yea 

Waahington.  D.  C. 

Yea 

Yea 

Yea 

Cleveland 

No 

No 

Yea 

Louiavllle 

Yea 

No 

Yes 

Yea 

New  Orleana 

Yea 

Yea 

Yea 

Burbank 

No 

No 

Yes 

Atlanta 

Yea 

Yea 

Yea 

Kanaaa  City 

No 

No 

Yea 

NOTE:  See  Note  for  Table  XXV 
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TABLE  XXVm 

COMPARISON  OF  PERFORMANCE  AND 
PATH  NOISE  SPECIFICATIONS  WITH  GLIDE  PATH 
FACILITIES  FOR  KC- 135,  DC-7,  TF-102 


95%  Prob 
Noise  Atnpl 
^  20  pa 


Facility 

MM  -Rwy 

KC-135 

DC -7 

TF-102 

Conflict 

Atlantic  City 

No 

- 

- 

- 

St.  Louis 

No 

Yes 

Yes 

Yes 

Yes 

Burbank 

No 

No 

No 

No 

San  Fran. 

Yes 

Yes 

Yes 

Yes 

Ontario 

No 

No 

No 

No 

NOTE:  See  Note  for  Table  XXV 
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Figure  46  Recommended  Course/Path  Noise  Amplitude  Criteria 
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3,  2.  3  Use  of  the  Recommended  Course/Path  Noise  Amplitude  Criteria 

As  a  result  of  the  preceding  studies,  Figure  46  shows  the  relation¬ 
ship  between  the  maximum  (96%  probability)  course/path  noise  amplitude 
and  the  distance  from  the  runway  threshold  that  is  being  recommended  in 
this  report.  If  the  noise  is  to  be  evaluated  in  any  region  of  the  approach 
the  theodolite-corrected  flight  recording  should  be  analyzed  for  a  time  in¬ 
terval  not  less  than  +10  seconds  and  not  more  than  20  seconds  about  the 
midpoint  of  the  region  to  be  evaluated. 

The  95%  maximum  amplitude  specification  is  the  allowable  percent¬ 
age  of  the  total  time  interval  in  which  the  course/path  noise  amplitude  must 
be  less  than  the  amount  specified  in  Figure  46  for  the  region  being  evalu¬ 
ated.  Figure  47  presents  a  typical  example  of  the  method  that  can  be  em¬ 
ployed  to  evaluate  the  course/path  noise  at  a  particular  facility.  If  the  sum 
of  the  time  intervals  t. ,  that  the  given  specification  is  exceeded,  is 

equal  to  or  less  than  5%  of  the  total  time,  f,  then  the  region  that  is  being 
evaluated  is  acceptable. 

Therefore.  100*'^  ~  *2  ^  “’^2  a:  95% 

T 


r 


0  -  THeoooure  coMpecreo  >sm^£  or  ^ouese/rv^rH  NOfse  -mo. 

0  -  mtoe  or  coutPse/frrM  Mf/se  ~m«- 
r  -  ffe^fON  ro  0e  er^tu^reo  (ttoroteo  secoNos  moe) 

±  Mjfsr  mrruruoe  sreariCRrtoN  rr  FotNr(^  -mo. 
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Figure  47  Recommended  Course/Path  Noise  Amplitude  Criteria 
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APPENDIX 

USE  OF  STATISTICAL  THEORY  IN  CLASSIFYING  ILS  FACILITIES 


1 .  SUMMARY 

This  appendix  discusses  the  classification  of  ILS  facilities  by 
statistical  decision  theory,  the  type  of  problems  the  theory  solves,  and 
the  validity  of  applying  an  optimization  criterion  of  the  theory  to  the 
classification  of  the  ILS  facilities  (acceptable  or  not  down  to  some  spe¬ 
cified  altitude). 

The  motivation  for  considering  the  applicability  of  decision  the¬ 
oretic  criteria  to  facility  classification  stems  from  the  initial  choice  of 
the  beam  error  power  spectrum  as  the  measurements  on  which  the  class¬ 
ification  of  each  facility  is  to  be  based.  In  effect,  the  choice  of  the  beam 
error  spectrum,  assigns  to  it  the  role  of  an  observable,  in  the  context  of 
decision  theory,  which  is  to  be  used  to  predict  the  conditions  of  all  the  land¬ 
ing  parameters  by  which  the  descent  operation  is  defined  as  acceptable  or 
not.  The  statistical  character  of  the  beam  error  spectra,  resulting  from 
measurements  for  each  of  25  ILS  facilities  presently  available,  was  found 
to  be  vinsatisfactory,  that  is,  the  spectra  turned  out  to  be  uncorrelated  with 
the  values  of  the  landing  parameters  as  found  by  exact  simulation  of  the  air¬ 
craft  control  system  fed  by  the  error  histories  of  each  of  the  facilities.  Al¬ 
though  it  can  easily  be  demonstrated  theoretically,  that  a  zero  beam  error 
spectrum  results  in  zero  landing  errors,  it  can  also  be  shown  that  the  beam 
error  spectral  measurements  are  uncorrelated,  nontrivially,  with  the  land¬ 
ing  errors. 

Following  the  demonstrated  failure  of  the  beam  error  spectrum  to 
correlate  with  the  simulated  landing  errors,  a  general  examination  was  un¬ 
dertaken  of  the  applicability  of  statistical  decision  methods  to  the  problem 
of  facility  classification.  The  conclusion  is  that  the  only  chance  of  obtain¬ 
ing  satisfactory  correlation  between  an  observable  and  the  true  landing 
parameter  is  provided  by  a  simulation  channel  approximation  of  each  of  the 
landing  parameters  for  its  corresponding  observable,  thus  giving  up  the  no¬ 
tion  that  frequency  distributed  measurements  on  the  input  parameter,  i.  e.  , 
the  beam  error,  can  correlate  with  each  of  the  landing  parameters.  This 
one-to-one  correspondence  between  observable  and  landing  parameter,  ad¬ 
mittedly  requires  more  complex  measurement  equipment,  and  one  can  safely 
assume  a  trade-off  between  the  order  of  approximation  (and  thus  the  com¬ 
plexity)  of  the  measurement  channels,  and  the  degree  of  correlation. 
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2. 


ELEMENTS  OF  STATISTICAL  DECISION  THEORY 


This  topic  is  discussed  following  consideration  of  the  type  of  con¬ 
ditions  under  which  an  ILS  facility  is  to  be  classified  as  acceptable  or 
rejectable  by  the  FAA.  A  sound  method  for  this  binary  classification  of 
future  facilities  is,  for  the  purposes  of  this  report,  the  objective  of  the 
study.  A  summary  of  the  conditions  and  their  relation  to  statistical  de¬ 
cision  theory  follows. 

First,  the  FAA  imposes  a  specification  on  several  descriptive 
landing  parameters,  such  as  errors  in  longitudinal  and  lateral  displace¬ 
ments,  altitude  rate,  pitch  and  roll,  these  errors  being  specified  at 
touchdown,  or  at  any  altitude  at  which  a  transition  from  automatic  to 
manual  landing  is  of  interest.  In  addition  to  the  aforesaid  instantaneous 
parameters,  a  mean  square  time  average  specification  imposed  on  control 
surface  activity  (mean  square  aileron  elevator  and  rudder  deflections 
over  the  duration  of  the  ILS  descent  operation)  would  be  consonant  with 
what  pilots  are  tactually  made  aware  of.  The  specifications  on  the  bounds 
of  these  errors  are  dictated  by  experience,  and  are  given  conditions  of  the 
problem  from  the  point  of  view  of  statistical  decision  theory. 

The  mathematical  theory  enters  into  the  problem  of  classifying 
an  ILS  facility,  when  the  information  available  to  tha  decision  maker  is  no 
more  than  statistical,  and  the  decision  maker  wants  to  find  decision  rules 
such  as  to  minimize,  or  bound  in  the  long  run,  his  misclassification  er^ 
rots. 


An  example,  taken  from  Reference  11  best  illustrates  the  problem: 

Suppose  a  weatherman  is  to  decide  each  day  between  two  hypotheses: 
Hq  -  it  will  rain  tomorrow,  and  Hj  -  it  will  be  fair  tomorrow.  The  decision 
is  to  be  made  on  the  basis  of  a  single  observed  quantity,  the  average  rate  of 
change,  x,  of  the  barometric  pressure  during  the  past  24  hours.  From  rec- 
rods  compiled  over  years,  the  probability  distribution  density  (p.  d.  d.  ) , 
p(x|Hg) ,  describing  the  distribution  of  x  on  days  preceding  rain,  and  the 


p.  d.  d.  ,  p(xl  H^)  ,  describing  the  distribution  of  x  on  days  preceding 
fair  ones,  are  p'recofnputed 

A  plausible  realization  of  these  two  distributions  is  illustrated  in  Fig¬ 
ure  4;-i  ,  with  the  24  hour  average  barometric  rate  on  days  preceding 
rain  being  distributed  more  negatively  than  on  days  preceding  fair  weather 
days. 


The  weatherman  can  now  utilize  this  statistical  information  by 
adopting  a  decision  rule  of  the  following  type: 

(a)  Predict  a  rainy  day,  if  the  observed  average  rate  x  is  less 
than  some  threshold,  x,  decided  upon  in  advance. 

(b)  Predict  a  fair  day,  if  the  observed  rate  x  is  greater  than  x. 

Now,  owing  to  the  overlap  of  the  distributions  j^x  IHq)  and  p(x  jHj) ,  no 
matter  what  choice  for  the  decision  threshold,  x,  mat  the  weatherman 
adopts,  he  risks  making  each  of  two  kinds  of  misclassifications; 

(c  )  Predicting  fair  weather  when  in  fact  it  will  rain, 

(d)  Predicting  rainy  weather,  when  in  fact  it  turns  out  to  be  fair. 


(In  general,  if  a  statistical  predictor  is  called  upon  to  predict  n 
categories,  rather  than  just  two  as  in  this  example,  he  risks  making 
n(n-l)  kinds  of  misclassification.  ) 


For  a  choice  of  some  particular  x,  as  a  point  which  divides  the  ob¬ 
servation  parameter  into  the  two  regions,  the  probability  of  making  the 
mistake  stated  in  (c)  is  then 


(1)  P 


H^ydx,  =  Qo(x) 


and  for  the  mistake  stated  in  (d) 


(2) 


Hjldx 


Q]  (x) 
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2.  1  BAYES  CRITERION 


Now,  if  the  weatherman  has,  in  addition  to  the  statistical  information 
described  by  the  conditional  p.  d.  d. ' s,  p(x|Hq)  andp(x|H2)i  (1)  the  prior 
probabilities  of  the  occurrence  of  rainy  and  fair  days;  that  is  independent 
of  the  observations  of  barometric  pressure  rates,  he  notes  the  relative  fre¬ 
quency  of  occurrence  of  rainy  and  fair  days,  and,  where  this  relative  fre¬ 
quency  of  occurrence  is  fairly  constant,  assigns  prior  probabilities 

(3)  P  {Ho}  =  W,  p(hi]=  1-p(Ho}  = 

and  (2) ,  a  realistic  assessment  of  costs  or  risks  on  some  scale  for  erring 
in  each  of  the  two  ways  described,  C^q  and  Cq^,  he  can  adjust  x,  the  parameter 
which  specifies  the  decision  strategy  in  this  case,  such  as  to  minimize  his 
average  risk, 

(4)  p[hi|h^}  p{h4cio  4  p(h„|hi}  Coi  =  c 

The  optimum  adjustment  of  x,  following  from  the  required  informa¬ 
tion,  is  put  in  evidence  by  noting  that  P^H2|  H  |  and  P  Hjj  are 
functions  of  x,  as  given  by  (1)  and  (2),  and  tha't  P  Hq  and  P  Hj_have,  by 
hypothesis,  known  values,  W,  and  1-W.  Hence  the  average  risk  C  is  a 
function  of  x,  and  the  cost  elements,  and  Cq^.  Equation  (4)  is  then 
rewritten 

(4a)  C(x,  Ro)  =  Qo(x)  W  +  Qj(x)  (l-WlC^i 


An  optimum  decision  threshold,  xg,  can  then  be  found  from  the 
constraint, 

(5)  ^  8et  =  0=  Q;,(xb)WCjq+  Qi'  {xb)(1-W)Coi 
which  yields 

(6)  XB  =  f(W,  Cqi,  Cjq). 
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For  the  decision  threshold,  x,  found  under  such  constraints,  the  weather¬ 
man  is  said  to  adopt  a  Bayes  strategy. 

2.  2  MINIMAX  CRITERION 

If  the  weatherman  does  not  know  the  prior  probabilities  of  occurrence 
of  rainy  and  fair  days,  W  and  1-W^  but  can  scale  the  costs,  C^q  and  Cq^,  and 
can  prescribe  the  two_conditional  d.  d. ' s,  p(xl  and  p{x  |H2)  ,  he  can  find 
a  decision  threshold,  Xj^,  which  minimizes  his  maximum  risk.  This  he  does 
by  considering  the  prior  probability,  W,  as  a  variable,  and  applying  the  con¬ 
straint, 

(7)  1=  0  =  Q,(VC,„-Q,(VCo, 

which  yields 

(8)  Xj^  =  giCjQ,  Cqj). 

It  may  be  noted  that  Xj^  depends  only  on  the  costs,  assumed  by 
hypothesis,  to  be  available. 

Now  if  the  prior  probability,  W,  is  considered  an  unknown  in  (4-a) , 
the  Bayes  c^straint  (5)  is  applied,  and  Xj^,  as  obtained  from  (7),  is  sub¬ 
stituted  for  xq,  in  (5),  the  result  will  be  a  constraint  on  the  unknown  prior 
probability,  W,  yielding  some  value,  i.  e.  , 

_ It  can  then  easily  be  shown  that  the  pair,  xj^,  is  a  saddle  point 

of  C(x,  W) ,  ^d  from  the  constraints  applied  to  obtain  and  it  can  be 

inferred  that  C(x]^,  is  a  minimum  with  respect  to  x,  and  a  maximum  with 

respect  to  W.  The  value  Xj^  as  obtained  from  (7) ,  thus  corresponds  to  a  Bayes 
strategy  under  the  most  unfavorable  prior  probability,  Wj^,  and  is  thus  optimal 
under  a  minimax  criterion. 

2.  3  NEYMAN  -  PEARSON  CRITERION 

If  the  weatherman  does  not  know  either  the  prior  probability  of  the 
occurrence  of  rainy  and  fair  days,  nor  the  costs,  C^q  and  Cq2  for  the  two 
kinds  of  error,  but  has  realistic  models  of  the  conditional  disturbances. 
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p(x'  Hjj)  and  p(x  H^),  he  obviously  cannot  minimize  his  average  risk, 

C,  in  any  sense.  He  can  however,  arrive  at  a  useful  optimization  of  x, 
the  decision  threshold,  by  imposing  a  constraint  which  practically  always 
has  significance  when  the  prior  probabilities  and  the  costs  are  not  avail¬ 
able.  This  constraint  consists  in  assigning  a  value  to  the  probability  of 
one  of  the  two  kinds  of  error,  i.  e.  , 

either  pJhj|hJ  =  a 

(10)  or 

and  minimizing  the  other  error  probability. 

For  the  example  chosen,  it  turns  out  that  no  minimization  is 
possible,  since  the  assignment  of  a  value,  a  tc^P|H||  t  or  ^  to 
P^HqI  determines  the  decision  threshold  x.  Another  example, 
however,  *of  wide  representation  will  be  shown,  where  both  an  imposed 
constraint  and  a  minimization  constraint  on  the  respective  errors  are 
required  to  adjust  x.  Before  taking  up  this  example,  a  few  remarks  on 
the  assigned  constraint  to  one  of  the  errors  is  in  order. 

Consider  first,  a  detection  radar,  which  is  required  to  make  auto¬ 
matically  the  binary  decision,  in  the  presence  of  interference;  target  or 

no  target  (H^) .  In  terms  of  the  output  of  the  detector,  this  amounts  to  de¬ 
ciding;  signal  plus  noise,  or  noise  alone.  Such  a  radar  would  require  as 
stored  information,  p(x|H2)  and  p(x  |  H^) ,  or  the  equivalent,  where  x  is 
now  the  output  voltage  of  the  detector,  and  p(x  |Hl)  and  p(x  |Ho)  the 
distributions  of  this  voltage  due  to  signal  plus  noise,  and  noise  alone.  The 
error  probabilities,  P^HqIHjI  and  pIHjIHq^,  correspond  now  to  the  pro¬ 
babilities  of  false  dismissal,  aftd  of  false  alarm,  respectively.  For  a  given 
environmental  signal  to  noise  ratio,  a  decrease  in  the  decision  or  alarm 
threshold,  x,  reduces  the  probability  of  false  dismissal,  but  must  concom- 
mitantly,  increase  the  probability  of  false  alarm,  as  may  be  seen  from  Figure  43, 
resulting  in  some  costly  action  being  taken  in  rain.  Under  such  conditions, 
the  policy  taken  is  to  impose  the  largest  false  alarm  probability,  or,  that  can 
be  tolerated,  arrived  at  from  considerations  external  to  the  statistical  cri¬ 
terion  discussed  here,  and  then  to  minimize  the  false  dismissal  probability, 
if  minimizable. 
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3.  THE  STATISTICAL  CLASSIFICATION  OF  ILS  FACILITIES 

(The  classification  method  discussed  in  this  paragraph  is  based 
on  the  Neyman- Pear  son  Criterion.  ) 

Assume  that  the  condition  of  a  single  landing  parameter,  y  (for 
example,  the  lateral  displacement  at  touchdown)  specifies  the  accept¬ 
ability,  Hq,  or  the  rejectability,  Hj,  of  the  ILS  facility.  More  specifically, 
in  terms  of  y. 

Ho.  if  -7  ^  y  ly^y 

Hp  ify<-y,  or  y  >  y  jyj>y" 


If  y  were  available  to  the  decision  maker,  there  would  be  no 
statistical  problem.  Owing,  however  to  the  fact  that  the  determination 
of  y,  for  a  given  facility  and  aircraft  control  system  would  require  a 
complex  aimulation  channel,  a  requirement  that  the  decision  maker  prefers 
to  avoid,  an  estimator  of  y,  is  determined  for  each  facility  and  aircraft, 
this  estimator  presumably  being  determined  from  a  substantially  simpli¬ 
fied  simulation  channel. 

The  decision  maker  now  is  confronted  with  prescribing  the  dis¬ 
tribution  densities: 

(11)  p(y||yj^y)  ^  p(y'|Ho) 
and  *  _ 

p(7||y|  >  y  )  =  p(7  |  h^) 


As  to  how  these  two  functions  can  be  found,  is  admittedly  a  basic  problem. 
First,  the  distributions  concern  a  future  population  of  facilities  for  which 
only 9^,  not  the  exact  parameter,  y,  will  be  available.  Let  this  population 
consist  of  N  members. 

Second,  a  pilot  population  consisting  of  M  members,  none  of  which 
is  in  the  N  member  population,  is  subject  to  the  dual  determinations,  y^^ 
and  yj^,  k  =  1  ,  M,  this  pilot  population  being  presently  in  existence,  and 

for  which  both  complete  channel  simulation  determinations,  y|^,  are  feasible. 
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Figure  49  Empirical  Correlation  Graph 
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The  two-dimensional  points  (yj^,  can  be  plotted  on  a  graph 
representing  the  y-y  plane.  (A  fictitious  example  is  shown  in  Figure 
49  ) .  The  estimator  axis,  'y,  can  then  be  partitioned  by  a  convenient 
uniform  class  interval,  The  collection  of  points  are  now  divided 

into  two  classes:  those  with  y  coordinates  inside  the  interval  ( -y ,  y) 
corresponding  to  the  H^-class,  and  those  outside  this  interval,  corres¬ 
ponding  to  the  Hj  -class.  Abstracted  then  from  the  M  member  population 
are  two  histograms: 


and 


(12)  h(y 


h(y 


II 


<:  y  ) 


y|>y  )  = 


mo(yi) 


MoAy 


(yi) 


MjAy 


m 


oi 


M  Ay 
o 


mii 


MjAy' 


where  moCy'^)  are  the  number  of  points  in  the  i*^  class  interval,  con¬ 
ditional  ujwn  their  y-coordinates  being  inside  the  acceptance  zone,  i.  e.  , 
for  |y  1j|  ^  y^  m^  (y^)  are  the  nvunber  in  the  i*-^  class  interval  conditional 
upon  •jyiil>y,  and  and  Mj  are  the  number  of  points  inside  and  outside 
the  acceptance  zone,  respectively,  and  necessarily. 


(13)  M  =  Mq  +  Mj. 


Deferring  temporarily  considerations  as  to  how  well  the  em¬ 
pirically  determined  histograms.  ( 12) ,  are  representative  of  the  fu¬ 
ture  population  distribution  as  required  by  (1 1) ,  it  will  now  be  shown 
how  the  Neyman-Pearson  Criterion  would  be  applied,  if  the  pilot  histo¬ 
grams  were  representative. 


81 


We  first  recognize  that  our  problem  now  is  that  of  determining 
two  decision  thresholds,  and'72  ob8_ervable,  y,  which  in  effect 

will  map  a  true  finite  acceptability  region,  _y  4,y<  y,  onto  a  finite  ac¬ 
ceptability  region  in  the  observation  space,  ^  shown  in 

Figure  50  >  the  true  classification  bounds  (— y,  y  ),  and  the  decision 
bounds  ,  divide  the  y-^  plane  into  four  regions  corresponding 

to  the  four  possible  outcomes  in  a  binary  hypothesis  test. 

We  must  also  decide  whether  to  fix  the  false  rejection  probability 
and  minimize  the  false  acceptance  probability,  or  the  other  way  around. 
For  tie  application  of  this  criterion  to  radar,  the  analog  of  the  former 
constraint  (fixed  false  alarm  probability-minimized  false  dismissal  pro¬ 
bability)  is  always  applied,  following  from  the  prevailing  operating  con¬ 
dition  being  that  a  detection  radar  is  overwhelmingly  not  presented  with 
any  targets. 


For  the  landing  system  situation,  we  must  expect  a  high  dominance 
of  safe  or  acceptable  landings  to  unacceptable  landings  (corresponding  to 
no  target,  and  target,  respectively  for  the  detection  radar),  otherwise 
we  have  a  landing  system  not  worth  taking  seriously.  Accordingly,  we  fix 
the  false  rejection  probability  and  minimize  the  false  acceptance  probab¬ 
ility.  For  the  observation  acceptance  zone  (^j  •^7*^92)’  utilizing  the 
empirically  determined  histograms,  there  results  the  false  rejection  pro¬ 
bability,  ^ 


(15) 


)dy 


^  h(y||y($y)dy 
Tz 


^  Qo(7j.T2)  =  'lo^M'  ^2)  ® 


(where  a  is  an  upper  bound  tolerable  to  the  decision  maker) ,  and  the 
false  acceptance  probability, 
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Figure  50  Conditional  Distribution  Densities 
Double  Decision  Threshold 
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1 


(16) 


>  y 


)dy 


m. 


li 


A 

=  Qi(yi.  y'2)  =  qi(ij.  i^) 


Using  the  method  of  l^gr^ge  multipliers  for  finding  the  pair  of 
optimal  decision  thresholds  (yj,  y2)  which  will  satisfy  the  value,  a, 
in  (15) ,  and  also  minimize  (16),  the  function  U  is  formed. 


(17)  U(yi,  72)  ^ 
conditions  are  imposed, 


^Qo(yi*  yz) »  and  the  following 


^U(yi,y2) 

3^1 


set  =  0 


aQi 

ari 


(yi.y2)  + 


(18) 


au(ri.y;) 
a  yz 


set 


0  ^ 


aoi 

dvz 


(y^.Tz) 


and,  in  addition 

Qo(^’  T2)  •  as  given  by  (15). 

The  first  two  equations  of  ( 1  8)  can  be  found  in  terms  of  the  given 
histograms  as 


(18a)  0=  -h(f;||y|>y  )  +  KH(y;|jy^^y  ) 

0=  h(y^Hy|>7  )  -\h(y'2|M^y  )• 

and  from  (IZ)  leads  to 
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(18b)  0  = 


mi 


m. 


M^Ay" 


0  = 


mii 


MjAy 


-X 


mo 


L2_ 


M^Ay- 


where  m. 


m. 


are  the  number  of  samples  in  the  acceptable  distribu¬ 


tion  in  claVs  inte^A^als  containing  the  decision  thresholds,  "^i  and'^2> 

>^ill  and  mi-  are  the  corresponding  numbers  for  the  unacceptable  dis¬ 
tributions.  i^om  (15)  and  (18b),  the  unknowns  X,  ii,  and  i2i  can  be  solved, 
which,  in  combination  with  the  histogram  tables,  ( 12) ,  leads  to  the  optimal 
decision  thresholds,  "y^  and^2- 

The  generalization  to  a  multiple  measurement  prescription  for 
acceptability  and  nonacceptability  is  straightforward. 

Let  yi ,  y2, ....  y  ,  be  n  measurement  parame^rs  for  which  ac¬ 


ceptability  is  prescribed*\»y  the  compound  event  jyij 
j/nl  ^  Yn'—****^  rejectability  by  the  compound  event 


^  yiand,  . 
yj  >71.  ot. 


and 


With  yj . .  y^j  as  tne  estimators  ot  yi ,  .  .  .  ,  y^  v 

obtained  distributions  in  a  pilot  population  can  be  obtained: 


^  as  the  estimators  of  yi , 


two  empirically 


^<>^1 . I  |yi(^  Yl .  jy.^^  Yn) 

^(Yl . Tn  I  |yi|  >  Yl  >  o*'  •  •  O'  |yn|  >  Yn^  * 

a  constraint  is  applied  to  the  false  rejection  probability. 
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.yii 

^  yin 

L 

dri...|drn  h(yp. 

r 

•••yn 

yp- 

•  |yn|^yn) 

j 

'y2i 

•  J  d9;h(9'p... 

y2n 

••^n 

yp' 

■  yn|  ^  ^n) 

and  the  false  acceptance  probability 

VZl  yZn 

Qi^ni^yzi . y'lw'^zn)  ••/  . yJi 

^11  yin 

is  minimized  via  2n-l  Lagrance  multipliers  in  2n  minimization  constraints, 
leading  to  2n  optimal  decision  thresholds,  p  ^\2‘  •  ’^^nl  ’  ^n2' 
measurements. 


|Vl|. 


>yi 


or 


'n 


>yr 


4.  THE  UNRELIABILITY  OF  CLASSIFYING  THE  ILS  FACILITIES  BY 
STATISTICAL  DECISION  METHODS 

As  stated  previously  in  this  report,  a  careful  investigation  into 
the  possible  correlation  between  power  spectral  density  measurements 
in  the  geometric  beam  error  and  any  of  the  landing  parameters  over  a 
25  facility  ensemble  (this  ensemble  being  subjected  to  complete  simula¬ 
tion)  failed  to  reveal  any  significant  correlation.  Specifically,  in  order 
to  reduce  or  eliminate  the  very  serious  error  of  false  acceptance,  one  is 
led  to  apply  decision  thresholds  to  the  spectral  measurements  (here  given 
the  role  of  estimator)  so  narrow  as  to  result  in  an  intolerably  high  false 
rejection  probability.  Note  that  this  is  a  result  possible  from  an  estimat¬ 
or  which  indeed  leads  to  zero  landing  error  when  the  estimator  is  zero, 
but  which  otherwise  is  not  correlated  with  the  parameter  for  which  it  has 
been  adopted  as  an  estimator.  It  is  not  difficult  to  see  why  the  course 
noise  power  spectrum  is  so  trivially  correlated  with  the  landing  errors. 

The  power  spectrum  (or  spectral  density),  for  any  spectral  fre¬ 
quency,  is  a  mean  square  over  a  time  interval.  But  a  bounding  of  a  time 
average  mean  square  really  leaves  the  course/path  noise  function  unspe¬ 
cified  except  in  the  trivial  case  of  imposing  a  zero  mean  square,  in  which 
case  the  noise  function  must  be  zero  everywhere  in  the  interval.  Any  land¬ 
ing  parameter,  on  the  other  hand,  must  depend  on  the  entire  history  of  the 
course/path  noise  preceding  the  time  instant  at  which  the  landing  parameter 
is  of  interest,  the  late  history  being  more  heavily  weighted  than  the  early 
history.  It  would  be  very  surprising  if  such  different  operations  acting  on 
a  given  ensemble  of  functions  were  statistically  correlated. 

From  these  considerations,  it  is  natural  to  seek  estimators  of  the 
landing  parameter  which  take  cognizance  of  the  simulation  channel  as  a 
point  of  departure,  and  which  then  approximate  the  weighting  characteris¬ 
tic  or  dynamics  of  the  landing  parameter  channel  by  means  of  a  significantly 
simplified  version  of  the  exact  simulation  channel.  This  approach  was  tried 
on  the  lateral  deviation  channel,  the  exact  channel  being  17th  order  and  slightly 
nonlinear,  the  approximate  channel  being  a  linear  order  filter.  The  re¬ 
sults  showed  more  correlation  than  for  the  power  spectrum,  but  nevertheless 
too  weak  to  be  used  for  statistical  prediction. 
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This  result  leads  naturally  to  considering  higher  order  approximations 
of  the  simulation  channels.  In  any  case,  this  point  of  view  requires  that  each 
landing  parameter  be  estimated  by  means  of  its  channel  approximation,  thus 
giving  up  the  notion  that  any  single  processing  operation  on  the  input  course/path 
noise  can  result  in  adequate  correlation  between  its  output  and  all  the  landing 
parameters. 

Accepting  the  notion  of  a  channel-wise  estimation,  we  arrive  at  a 
fundamental  problem  concerning  the  statistical  correlation  of  an  estimator 
with  the  parameter  it  is  to  estimate,  and  the  confidence  that  can  be  put  in  em¬ 
pirical  distributions  derived  from  a  small  sample  population. 

If,  in  the  absurd  limit,  the  estimation  channel  were  an  exact  replicate 
of  the  simulation  channel  which  by  definition  determines  the  true  landing 
parameter,  no  trials  would  be  necessary  for  obtaining  the  two  conditional 
distributions,  these  distributions  being  required  for  attaining  Neyman-Pearson 
adjustments  of  the  estimator  decision  thresholds.  The  estimator  decision  thres¬ 
holds,  in  this  limiting  case,  would  be  exactly  the  FAA  specified  boiinds  on  the 
acceptance  zone  of  the  parameters  in  question. 

If  now,  at  the  expense  of  high  order  approximation  for  each  estimation 
channel,  a  high  degree  of  correlation  is  attained  in  the  pilot  experiment,  it 
can  be  concluded  that  the  future  population  will  also  be  highly  correlated,  and 
that  the  two  distributions  empirically  obtained  from  the  pilot  population  can 
serve  with  confidence  as  the  corresponding  distributions  for  the  future  popu¬ 
lations.  Stated  in  other  terms,  a  strongly  correlated  estimator,  results  in 
only  slight  overlap  of  the  two  distributions,  a  characteristic  which  is  more 
likely  to  be  extrapolable  from  a  small  to  large  population  than  for  a  large  over¬ 
lap.  Hence,  there  arises  the  question:  Given  the  two  distributions,  empirically 
derived  from  a  small  population,  and  not  necessarily  with  small  overlap,  are 
the  probabilities  of  misclassification  (false  rejection  and  false  acceptance)  re¬ 
sulting  from  a  Neyman-Pearson  adjustment  of  the  decision  thresholds  of  the  es¬ 
timator  for  the  pilot  population  necessarily  close  to  or  not  close  to  the  corres¬ 
ponding  misclassification  probabilities  for  the  large  future  population?  This 
question  has  not  been  resolved  at  present.  However  it  seems  very  doubtfixl 
that  small  population  estimation  distributions  which  (1)  result  from  poorly  corre¬ 
lated  estimators,  and  which  (2)  are,  for  the  ILS  facility  classification  problem, 
multivariate  distributions,  will  be  necessarily  close  to  the  large  and  decisive 
problem  distributions. 

Moreover,  since  it  may  be  the  intention  that  only  these  25  facilities 
or  small  sampling  of  the  total  number  of  facilities  will  be  used  for  lower 
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altitude  approaches,  there  is  no  point  in  using  statistical  decision  methods 
for  classifying  these  facilities.  The  implications  of  this  statement  are 
either  (1)  investigate  higher  and  higher  order  estimation  channels  until  very 
strong  correlation  is  found  between  the  outputs  of  the  estimation  and  exact 
channels  or  (2)  classify  the  activity  on  the  basis  of  exact  simulation  only. 

In  view  of  the  complications  of  Item  (1)  with  its  nonelimination  of  misclass- 
ification  error,  Item  (2)  is  recommended  as  the  technically  best  method  to 
be  followed. 
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